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1. Introduction 


This report describes an analysis of the electromagnetic (EM) fields which are 
backscattered from a two dimensional S-shaped inlet cavity. Figure 1.1 shows the 
geometry of this problem along with all the variable parameters. The interior 
walls of the S-shaped cavity are perfectly conducting with a thin absorber coating. 
The planar interior termination is made perfectly conducting. The S-shaped inlet 
is assumed here to be modeled by joining appropriate sections of parallel plate 
waveguides and uniform annular waveguides as is evident from Figure 1.1. 

The two basic scattering mechanisms of the inlet cavity of Figure 1.1 are the 
scattering from just the edges at the open end, and the scattering from the interior 
termination within the inlet cavity. A plane wave is assumed to be incident on 
the inlet at some aspect angle (0) and encounters the edges at the open end of the 
inlet. The two leading edges thus scatter part of this incident energy by diffraction at 
those edges to give the direct contribution from the opening by itself. The scattering 
from the interior cavity region is due to the incident energy which couples into the 
natural modes of the piecewise uniform waveguide sections of the S-shaped inlet, 
and propagates to the end termination to be then reflected back into the aperture 
at the open end from which it radiates into the exterior region. In addition to 
these primary scattering mechanisms, there are higher order effects which need 
to be considered. These include the small reflections due to the discontinuities at 
junctions between sections, the slightly larger internal reflection of the energy by the 
aperture, and the multiple wave interactions between all these scattering centers. 

Chapters 2 and 3 describe an efficient hybrid combination of asymptotic high 
frequency (HF) methods with the modal technique that is used to analyze the 
scattering from a perfectly conducting inlet [1,10]. The fields in the guide are 
expanded into the natural waveguide modes of each uniform section. The relative 
strengths of these modes are represented with coefficients which can be expressed 
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in column matrix form. Reflection at a junction, transmission through a junction 
and propagation through a section are all represented as rectangular matrices which 
operate on the column matrices for the coefficients of the modes incident on the 
junction. The elements of these matrices are found using high frequency techniques 
such as the Kirchhoff approximation, the equivalent current approach, and modal 
ray-optics. The multiple interactions are accounted for using the self-consistent 
Multiple Scattering Matrix (MSM) formulation. The latter is essentially the same 
as the Generalized Scattering Matrix formulation [2]. 

Chapter 4 extends the method of chapters 2 and 3 to include inlets with thin 
absorber coatings on the inner walls using a small perturbation to the propagation 
constants of the waveguide sections. This method is limited to relatively small 
values of absorber thickness and loss. 

Chapter 5 describes a numerical method of finding the modal reflection coef- 
ficient for the reflection from the first junction of the inlet geometry. Refering to 
Figure 1.1, this is the reflection of a waveguide mode of section 1 incident at the 
junction between sections 1 and 2. The method expands the fields in each waveg- 
uide section into modes and solves for the unknown coefficients by matching the 
fields at the junctions. It is not very efficient but it can handle arbitrary impedance 
boundary conditions on the waveguide walls. This solution is included to test the 
modal reflection obtained from the perturbational approach of Chapter 4. 

Chapter 6 describes a purely ray-optical method of analyzing cavity scattering. 
The geometrical optics (GO) incident field which enters the inlet from the aperture is 
then tracked via rays through the inlet as they undergo multiple reflections from the 
inlet walls to reach the termination and then bounce back again to the aperture. An 
aperture integration (AI) is subsequently performed on this GO field scattered back 
into the aperture to obtain the exterior far field radiation from the aperture. The 
advantage of the GO/AI method is that it can handle arbitrary inlet geometries 
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other than inlets made up of uniform waveguide sections. Also, it can handle 
any type of absorbing wall material which can be characterized locally by plane 
wave reflection coefficients. However, this approach which ignores diffracted rays 
becomes less accurate for electrically small inlets and/or long inlets. It is included 
in this report mainly for comparison purposes and to supply an alternative analysis 
method. Also, this GO/AI technique, which ignores diffraction effects at the open 
end that also propagate into the inlet, is found to provide a reasonable estimate of 
the average of the backscattered field but not the details in the field pattern. 

A few numerical examples along with some measurements are included in chap- 
ters 4, 5 and 6. Throughout the report an e JU,t convention is assumed and suppressed. 
“TE” (Transverse Electric) means that the E-field is transverse to the direction of 
propagation, or equivalently, the E-field has only a component perpendicular to 
the plane of incidence. Similarly, “TM” (Transverse Magnetic) means the H - field 
is perpendicular to the plane of incidence. Here, the plane of incidence is defined 
by the incident ray and the axis of the first parallel plate waveguide section of the 
S-shaped inlet. 
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2. Self-Consistent Multiple Scattering Matrix Formulation 


The inlet considered in this report is made up of interconnected piecewise sec- 
tions of parallel plate waveguides and circularly curved or “annular” waveguides. 
The fields inside these guides can be expressed as an expansion of only the prop- 
agating ortho-normal waveguide modes. 1 Ortho-normal simply means that each 
mode is normalized to carry unit power and carries it independently of all the other 
modes in the guide. The modes themselves satisfy the wave equation and boundary 
conditions of the corresponding uniform guide which extends to infinity in both 
directions. For the annular guide, this would mean an infinitely extended angular 
space. To account for the finite lengths of the various guide sections and the multi- 
ple interactions between all the junctions and the open end, the Multiple Scattering 
Matrix (MSM) formulation is used [2]. 

A junction between two dissimilar waveguide sections can be characterized in 
terms of transmission and reflection matrices associated with the junction which 
relate the coefficients of the waveguide modes on either side of the junction. Figure 
2.1 shows a typical junction between waveguide sections which are both semi-infinite 
(in this case, it is the junction between sections 1 and 2). Energy propagates down 
section 1 from the left and is incident on the junction. In general, a junction is 
formed by joining a p th waveguide section with a q th waveguide section. Thus, 
the sections p and q are physically connected together at the junction. Using the 
symbol U for the ^-directed JE7-field of the TE case and the //-field of the TM case, 
respectively, the incident, reflected and transmitted fields at any junction between 
the p th and q th waveguide sections can be expressed as 

U ine = £ A n u n e~^ nX ( 2 . 1 ) 

n=l 

1 Since the separation between the junctions is sufficiently large, the evanescent modes are ignored 
as they contribute negligibly. It is noted that evanescent modes cannot be normalized the same 
way as is done here for the propagating modes. 
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( 2 . 2 ) 


£?"" = jr B„u„e ie ” 

n= 1 
M 

jjtrant = £ C m U m j Vm * (2-3) 

m=l 

where u n represents the n th transverse ortho-normal modal field of section p at the 
junction and u m likewise represents the m th modal field of section q at the junction. 
The coefficients A n , B n , and C m can be represented as column vectors 


[A] 

[B] 


A a 

An 

Bx 

Bn 


(2.4) 


(2.5) 


[C\ = 



( 2 . 6 ) 


L Cm J 

The reflection and transmission scattering matrices for this junction between the 
p th and q th waveguide sections are then defined by 


[■ B } = 


[C] = [T '„}[A] 


[«„] = 

R'A 

p21 


pi N ■ 
•** K rx, 





[T„] = 

’ rpH 

P<2 

rp2\ 

rpX2 

T lN 1 

iSr 


rpMl 
l Apq 

rnM2 

rpMN 

Apq 


(2.7) 

( 2 . 8 ) 

(2.9) 


( 2 . 10 ) 


Here, the reflection matrix is associated with the fields reflected back into the p 1 


th 


section when it is joined directly to the q th section. Likewise, the transmission 
matrix is associated with the fields transmitted into the q th section when the fields 
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are incident from the adjoining p th section. N and M are the number of propagating 
modes in the adjoining sections p and q , respectively. Notice that [itp,] is N x N 
and [Tpg] is M x TV. 

For the special case of the open end (the junction between waveguide section 1 
and the exterior region) the scattering matrices have a slightly different form. For 
the coupling into the inlet and the scattering from the open end, the fields can be 
expressed similarly to equations (2.1), (2.2) and (2.3) as 

U inc = ^U ince -M^co.f+ynm0) ( 2 . 11 ) 

U rtfl = ( 2 . 12 ) 

y/p 

N 

fjtran. _ C n U n e~ ik ^ X (2.13) 

n=l 

where U inc is a scalar representing the magnitude of the incident plane wave and 
p is the distance to the receiver referenced from the lower (*, y = 0) edge. U Te * 1 
represents the field scattered back into the exterior by the edges at the open end, 
and U trani represents the field coupled (or transmitted) into the inlet. B is a scalar 
for this case and [ C ] is a column matrix as before. The reflection and transmission 
matrices can now be defined similarly to equations (2.7) and (2.8) as 

B = Ro!U ine (2.14) 

[C] = (2.15) 

where the subscript “0” refers to the exterior region. Notice that the reflection 
matrix has reduced to a scalar for this case and the transmission matrix has reduced 
to a column matrix (because [C] is a column matrix). The matrices describing the 
exterior radiation and the interior reflection from the open end when modal fields 
are incident from inside the guide are defined similarly as follows. 

u inc = Y, A nU n e ik0nX (2.16) 

n—1 
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Figure 2.2: The junction of Figure 2,1 with a termination in section 2. 


N 

U refl = B n^ne- jk0nX (2.17) 

n=l 

u tran ‘ = zC—-r (2.18) 

VP 

[B] = [* 10 p] (2.19) 

C = [r l0 ][4]. (2.20) 

Notice that [_R 10 ] is a square matrix and [T 10 ] is a row matrix (because C is a scalar). 

To include the effects of more than one scattering junction, the Multiple Scatter- 
ing Matrix (MSM) is used [2]. To illustrate its implimentation, consider the simple 
geometry of figure 2.1 with a termination placed in section 2 as shown in figure 2.2. 
The total reflection matrix [i?i 2 ] in region 1 due to the junction between waveguide 
sections p = 1, q = 2, as well as the termination is given via the MSM for that 
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junction as 




" sM i = [« 12 i + 

[r„l[ft][Jiil[ft]([i] - Pi.] 


j 


( 2 . 21 ) 


in which 

(jy = 


e ~ iu 0 

0 
0 


0 


( 2 . 22 ) 


0 e ~ iUM4, ° 

where [/] is an identity matrix, [P 2 ] describes the propagation through section 2, 
and [R t ] is a reflection matrix describing the termination. As an example, if the 
termination were a perfect electric conductor, for the TM case [fit] would equal [/] 
and for the TE case [il t ] would equal — [I]. The above result for [R^ SM \ accounts 
for all the multiple wave interactions between the junction and the termination. 

To describe the multiple interactions of more than 2 scattering junctions (such as 
in the inlet of Figure 1.1), the MSM result of (2.21) is modified by simply replacing 
[fit] with the MSM of the next junction [1]. This is repeated until the final junction 
associated with the termination is reached. For the inlet under consideration (see 
Figure 1.1) there are four junctions formed by connecting five dissimilar waveguide 
sections, a termination, and the open end which amount to a total of six different 
scattering junctions. Modifiying the MSM to include all six junctions, the total 
scattered field can be expressed as follows. 

,-ikp 


jjtcat jjinc^_ 


Vp 


r usm = fi 01 + [r 1 „HP 1 ](fiK s "i(P,]- 

• ([/] - [fiio]lPiPf$ ;5 *'|[Pi|)"' [r 01 ] 


(2.23) 


(2.24) 


= [flu] + Pm] [P s pg SM ][Pj] • 

• ([/]-[fl 2I ][P ! ][fl5| SM ][P,])" I [T 12 ] 


(2.25) 
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( 2 . 26 ) 


= [*«] + [r, 2 ](p 3 i(jj" SM i(p 3 ] ■ 

•(W-[fl3 J ][P3][fi5J s "||P s ))"'(r I 3] 

[*£“] = i* 34 i + • 

• ((/] - [fl<3](P,](^ SM !(/>,]) ■' |r„] (2.27) 

[^MSM] _ [^ 45 ] + [7’ s4 ]{P s ](i? ( |(A] • 

•((/]-[B,4|[P,)[«,)[P s ))- I (r«| (2.28) 

where the junction matrices are defined using the “pg” convention as before. 

It is seen from symmetry considerations that some of the above matrices are 
equal. Also, as will be shown later, the transmission matrices at a junction are 
reciprical. Appendix D lists these relations. With these simplifications, only eight 
of the original eighteen reflection and transmission matrices need to be evaluated. 
They are [R 12 ], [#21], [R34], [#43], [.Rio], [Rt], [2\s], and [T B4 ]. If the radii of curvature 
of the annular waveguide sections 2 and 4 are equal, i.e., if r 2 = r 4 then only [Ri 2 ] , 
[R 21 ], [Rio], [Rt], and [T 12 ] need to be evaluated. 
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3. Hybrid Asymptotic Modal Analysis for Determining the Elements 

of the Scattering Matrices 

3.1 Modal Field Structure 

To determine the elements of the scattering matrices, a knowledge of the modal 
fields of the two types of uniform waveguides must first be obtained. Figure 3.1 
shows a parallel plate guide which extends to infinity in both directions. The total 
field propagating in the ±i- direction inside the parallel plate guide can be expressed 
as a summation of the orthonormal waveguide modes of the guide: 


E £ {x,y) 


(3.1) 

H ± (x,y) 

= 'EAihiMe™-' 

n 

(3.2) 

fin 

= v**-<t>* 

(3.3) 


where the coefficient A* determines the relative strength of mode n. The orthonor- 
mal modes are normalized according to 

f o [cn(y) x h$(yj\ -(±x)dy = 1. 



Figure 3.1: Parallel plate waveguide geometry. 
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In any waveguide there are an infinite number of modes. However, there is 
always a finite number of propagating modes, the rest being evanescent (they decay 
exponentially along the axis of the guide). The evanescent modes are ignored in this 
analysis as they contribute negligibly for the reasons explained in chapter 2. If only 
the propagating modes are included, the summation is truncated at N — int(kd/ir), 
the integer portion of kd/ir. Notice that the number of propagating modes increases 
linearly with the width of the guide. Therefore, electrically large inlets may support 
too many waveguide modes to allow an efficient modal analysis. Chapter 6 discusses 
an alternative ray method useful for large guide widths, however, the latter is not as 
accurate because it ignores diffraction effects at the open end which propagate into 
the guide and can eventually dominate over the reflection effects over long paths. 

For the TE case 


e${y) = e n (y) 


= zP n sin(^y) 


hniv) = 


±M(y) + h no (y) 
PnY 0 


= / 


. ,nir « nir nir 


2k Zo 

d& 


(3.4) 

(3.5) 

(3.6) 


and for the TM case 


^n(y) = My) 


My) = 


P n = 


zP n cos(^y) 
±e„ t (y) + Cno(y) 
P n Z o 


. ,mr . . nir . .nir . 


! 2kYo 

dftnCn 


(3.7) 

(3.8) 

(3.9) 
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(3.10) 


£n 


2 if n = 0 
1 otherwise 


where Zo is free space impedance, Yo is free space admittance and k is the free 
space wave number, 2ir/\. The subscripts and “a” signify the transverse and 
axial components, respectively, of the fields which lie in the plane of incidence (in 
this case, the y and x components, respectively). Notice that for the TE case there 
is only a z-component of the E-field. Similarly, there is only a z-component of 
the iT-field for the TM case. This helps simplify the analysis of two dimensional 
structures. The derivation of these modes is given in Appendix A. 

These modal fields can be physically interpreted in terms of an equivalent set 
of modal ray fields by expressing the sine and cosine terms in their Euler form and 
combining their exponential arguments as shown in Appendix A. For the TE case 


£n{y)e 





jk(y sinflnT* c ° 8 *») 


e jfc( - V « n On =f * coa On ) 


(3.11) 


and for the TM case 


h n (y)e T ^ nX — z— e M - V «in «« T* C 03 0 n ) j 


(3.12) 


where 


0 n = iin-'Q. (3-13) 

As figure 3.2 shows, the above alternative expressions for e„(y)e ±J ^ nX and h n {y)e^^ nX 
represent oppositely traveling plane waves within the guide which make a charac- 
teristic modal ray angle of 0 n with the walls. This ray-optic interpretation will be 
important when the reflection from a junction is calculated later. 

Figure 3.3 shows an annular guide which extends to infinity in the indirections. 
This, of course, is physically impossible because it would just connect up with itself 
to form a closed ring. However, this is allowable mathematically, and is neccesary to 
avoid the periodic <j) boundary condition of an annular guide of 360 degrees because 
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the bent portions of the inlet are made up of sections of annular guides which are 
fractions of this. The total field propagating in the indirection can be expressed 
as (see Appendix B) 


EHpA) = XX 

n=l 
n= 1 

For the TE case 



hn(p) 

Rn(kp) 

A n 

R n (ka) 


Cn{p) 

zAnRnikp) 
±h nt (p) + hna ( p ) 
A n Y 0 


ip^Rnikp) - 4>jR' n (kp) 


Yj [Hll\kb)Hl'J{kp) - HH\kb)Hll\kpj\ 


a 


0 


(3-14) 

(3.15) 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


where the transcendental equation (3.20) is solved for the eigenvalues u n . For the 
TM case 


Kip) = Kip) 

= zAnRnikp) 

K(p) = ± e nt(/>) + ena(p) 

= A n Z 0 

K(kp) = i [v£>(W)H<l>(M - flS v (W)v£>(M] 


*i>fR„(kp) + <hK(kp) 


/U = 


du 

u 




( 3 . 21 ) 

(3.22) 

(3.23) 

(3.24) 
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R' n (ka) = 0 (3.25) 

where the transcendental equation (3.25) is solved for the eigenvalues u n . 

H^', and H^' are the Hankel functions of the first and second kinds of order u n 
and their derivatives with respect to argument, respectively. For this case p is the 
transverse direction and <f> is the axial direction. 

The modal ray form of the modal fields in the annular guide are obtained us- 
ing the Debye Tangent approximation for the Hankel functions. This is done in 
Appendix B. The geometrical interpretation is shown in figures 3.4 and 3.5. It 
consists of oppositely crossing cylindrical waves that form a circular modal caustic 
which is concentric with the walls of the guide with radius v n /k. The angles these 
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Figure 3.6: Field of a whispering gallery type mode, 
modal rays (with cylindrical spreading) make with the guide walls are given by 

= cos-'ljjjj) (3.26) 

0(m = COS _ '(^~). (3.27) 

As shown in figure 3.5, the caustic contour can lie inside the guide. This spe- 
cial case is called a “Whispering Gallery” (WG) mode and is encountered in most 
concave surface guided wave problems. This particular type of waveguide mode 
confines most of its power in between the caustic contour and the outer wall, as 
shown in figure 3.6. So, ray-optically, the propagating part of a WG mode never 
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“sees” the inner wall even though its eigenvalue is dependent on the existance of 
the inner wall. 

It is seen from the ray-optical forms that the waveguide modes of the parallel 
plate guide and the annular guide are closely analogous. In fact, the number of 
propagating modes are almost always the same for the two types of guides if they 
are of the same width. Furthermore, it is straightforward to show that the modes of 
the annular guide approach those of the parallel plate guide as the radius approaches 
infinity. 

3.2 Evaluation of the Elements of the Scattering Matrices 

The elements of all the scattering matrices are calculated by first replacing the 
discontinuities in the guide with equivalent currents which simulate the effects of 
the discontinuity [1]. The reflection and transmission type matrix elements are then 
simply the excitation coefficients due to these currents radiating in a uniform guide 
of infinite extent when a mode of unit amplitude is incident on the junction in 
question. For example, to find the elements of [ffo] the junction would be replaced 
by a uniform parallel plate guide with equivalent currents where the junction used to 
be (see Figures 3.7 and 3.8). It is then important to find these equivalent currents 
accurately. The excitation coefficient of the n th mode due to current sources in a 
uniform guide (see Figure 3.9) is given by [1] 

C n = [«« ’ J- Mnt + M * M\ ds (3.28) 

c i = (3.29) 

— * — ♦ 

for the TE and TM cases, respectively, in the ±a directions. J and M are the 
equivalent electric and magnetic current sources and a is the unit vector in the 
axial direction. The subscripts “f” and “a” again mean “the component in the 
transverse direction (to a)” and “the component in the axial direction”. It is noted 
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that (3.28) and (3.29) apply to both transmission and reflection calculations where 
the integrations are performed over the appropriate source regions which must lie 
in a cross section of the guide. 

3.2.1 Junction Transmission Matrices 

The transmission matrix elements T™ n are found by replacing the incident modal 
field with equivalent incident electric and magnetic surface currents in the aperture 
formed by the junction and then evaluating the excitation of the transmitted modes 
by these currents using (3.28) and (3.29). To illustrate the procedure, [T 12 ] will be 
derived for the TE case. [T 12 ] describes the transmission of energy from section 1 
across the junction to section 2. T™ is the transmission matrix element which is 
the excitation coefficient of the m th transmitted mode in section 2 when the n th 
mode of section 1 is incident on the junction. Using equation (3.28), it is given by 

T T = -\[ fcm(p) • /« - [h mt (p) + L(/»)| • Mn} dp (3.30) 

where a is <f> and 


e«(/>) 

= zA m Rm(kp ) 

(3.31) 

hmt(p) 

= pA m Y^R m (kp) 

(3.32) 

h ma (p ) 

= -ijA m Y 0 R' m (kp). 

(3.33) 


—4 — ♦ -4 *4 

The equivalent currents, J — n x H and M = E x n are given by 
J n = x x h+(y) 

= -zP n Toy sin(~ y) (3.34) 

M„ = e+(y) x x 

= yP n sin(^y). (3.35) 
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(3.36) 


Substituting back into equation (3.30) gives 

T S " = + 

Using y — p — b and changing variables to u = kp results in 

= \ a " p "t£ (j + t) 1 sin lv/ u - H < 3 - 37 > 

The integration is done numerically or possibly through the use of asymptotics if 
the guide is large enough. T$ n is found similarly for the TM case 

*3" = r (t + v) CM - H (3 - 38) 

If T™ is derived using the same procedure, it is found that T™ = T™, as men- 
tioned earlier in Chapter 2. T 5 ™ n is the same as T$ n with the appropriate changes 
in radii a and b. 

The equivalent sources J n and M n are based on the modes incident on the 
junction. This approximation should generally provide the dominant contribution 
to T™- A correction to this approximation can be deduced via a modification of 
Ufimtsev’s Physical Theory of Diffraction in which equivalent currents based on the 
fringe wave diffraction by the junction are used to obtain this correction [3]. The 
latter correction is usually weak and is therefore neglected here. 

3.2.2 Transmission Matrices of the Open End 

The coupling of the incident plane wave into the first section is described by 
the column matrix [Toi], as defined in chapter 2. The elements of [Toi] are found 
using, once again, equivalent currents and the excitation coefficient equation very 
similarly to the manner in which the junction transmission matrices were derived 
above. The element is the excitation coefficient of the n th mode of section 1 due 
to the plane wave of unit amplitude incident at angle 9. For the TE case 

T oi = + (3-39) 
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The equivalent currents J and M are found as before using the Kirchhoff approxi 
mation indicated previously in section 3.2.1, namely 


M = E inc x x 


ze~ jkvaiat x x 


~ e -jky sin* 


J = x x H inc 

= X x Y 0 e~ ikv ™ 9 (x sin 9 - y cos 9) 
= -zY 0 e- jkv,in * cos$ 
e n {y) = zP n sm(^jy) 

h n t(y) = -yPnYo^- sin(^y). 


(3.40) 


(3.41) 

(3.42) 

(3.43) 


Substituting back into equation (3.39) and integrating gives 
Toi = ^Y 0 P n de-rt kd ™ 9 j n - 1 ^cosO + (j^ 

• |sinc ^kd(sm9 n — sin0)j 
— (— l)"sinc ^fcd(sin #„ + sin*)] } 

where 


(3.44) 


sine* = . (3.45) 

X 

Notice that the sine function peaks when 9 = ±0„, which is understandable physi- 
cally when the modal ray form of the modes is considered. The coupling is maximum 
when the incident plane wave lies along a modal ray angle. 

The TM case is found similarly as 
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(3.46) 


+(— l) n sinc |^fcd(sin0 n + sin0)j | . 

The elements of the matrix [T 10 ] are found using the aperture integration (AI) 
technique. The element T™ 0 is the coefficient of the radiated field due to the n th 
mode of section 1 being incident on the open end. For the TE case, it is defined by 
e~i kp I ik e.~i kp r d r , 

[co s0Er{y) + ZoH?(yj\e-to™'dy (3.47) 

where p is the distance to the receiver referenced from the lower edge (x,y = 0) 
of the inlet. This integral is derived in Appendix E. The aperture field is the n th 
modal field, 


E? p (y ) = z-e n (y) 

= P n sin(^y) (3.48) 

H a v p {y ) = y-h~(y) 

= P„r 0 ysm(— y). 

Substituting for the aperture field, integrating and removing the cylindrical wave 
term gives the matrix element as 


rpn 

-MO 


1 /jk 

4V27T 


P n de- J ’ kd ' inf j n - 1 (cos0 + 


■ | sine y^fcd(sin 0 n — sin 0)J 
— ( — 1 ) n sinc ^ kd(sin B n + sin 0)j j . 


(3.49) 


Notice again that due to the sine functions, the radiation is maximum along the 
direction of the modal ray angles. Also notice the close similarity between the 
expressions for P 0 " and T” 0 . This is due to reciprocity. 

The TM case is found similarly and is given by 

r "« = 


I 



(3.50) 


* jsinc ^fcd(sin 0 n — sin 0)j 

+( — l) n sinc ^kd(sind n + sin0)j j . 


3.2.3 Reflection Matrix Elements for a Junction 

The reflection matrix elements are found by replacing the discontinuity of a 
junction with equivalent line currents placed on the walls at the junction and then 
finding the excitation of the reflected modes by these currents [1], The equivalent 
line currents are found using the Geometrical Theory of Diffraction (GTD). First, 
the modal ray form of the incident mode is used to find the modal ray field incident 
on the discontinuity. That incident modal ray field undergoes diffraction at the 
discontinuity. Then the line current which radiates in the presence of the waveguide 
walls is found using the diffraction coefficient for the discontinuity. The method will 
be illustrated by deriving [Rio] and [R 12 ]. 

The geometry relevant to [Rio] is shown in Figure 3.10 and the equivalent ge- 
ometry with the discontinuity replaced is shown in figure 3.11. R" o n is the reflection 
matrix element which is the excitation coefficient of the m th mode reflected from 
the open end of section 1 when the n th mode is incident from section 1 . For the TE 
case, using equation (3.28) 


nmn 

^10 — ~~ 

\J 0 {Ml/) ' ~ I ~M(y) + hma{y)\ ■ tin} 

(3.51) 

where a = x and 



My) = 

- n ■ , mlT i 

zPm sin(-^-y) 

(3.52) 

h mt (y) = 

-T, v Pm . ,rnir 
-yP m Y 0 ~- sm(~y) 

(3.53) 

hmaiy) ~ 

'■n v m7T / m7r i 

- xPmY °Wd co< ~d y) ' 

(3.54) 


Since the line sources are located at the junction on a conducting plane, the electric 
source is shorted out leaving only magnetic currents. For the TE case these magnetic 
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Figure 3.10: Geometry relevant for finding [i2 10 ]. 
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Figure 3.11: Aperture replaced by equivalent currents. 
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sources are a line of x-directed infinitesimal dipoles located at Q 0 and Qd- The 
equivalent sources are therefore 

J n = 0 (3.55) 

M„ = zM*(Q 9 )6{y) + xM?>(Q d )6(y - d) (3.56) 

where S(y) is the Dirac delta function. Substituting back into equation (3.51) gives 


i *r. n = 


2 ~ "‘~' 1 jkd 1 — n 1 v ~~ n ' (3.57) 

It remains only to find the scalars M* P (Q 0 ) and M* p (Q<j), which are the strengths 
of the line dipoles. This is done by finding the modal ray form of mode n incident 
on the discontinuities and using the appropriate diffraction coefficient [1] which is 
presented in Appendix C. The line source strengths are then given by [4] 




(3.58) 

(3.59) 


where 0 m and S n are the modal ray mode angles of the reflected and incident modes 
and, using equation (3.11), 


C(<?oM = ~j 

zrm-t = -(- 1 rfy. 


(3.60) 

(3.61) 


Substituting equations (3.58), (3.59), (3.60) and (3.61) back into equation (3.57) 
gives the final expression 

-P m P n Y Oy J^D,{0 m ,9 n ) iim + n even ^ ^ 

0 otherwise 

[i? 10 ] for the TM case is found similarly, except that magnetic line sources are 

used rather than magnetic line dipole sources [4], 


R? n 


no 


■{ 


M n = zM“ n (Q 0 )6(y) + zM“ n (Q d )6(y - d) 


(3.63) 
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Figure 3.12: Geometry relevant for finding [/Z 12 ]. 


MH"(Q d ) = -Zo^D H (« m ,e n )h?(Q d )-z. 


The final expression is then 


R 


mn 

10 



m + n even 
otherwise 


(3.64) 

(3.65) 


(3.66) 


Notice that since the diffraction coefficient is reciprical, [i? 10 ] is symmetric. This 
reduces the number of computations considerably. 

Figure 3.12 shows the relevant geometry for finding [i2i 2 ]. The matrix [iZ 12 ] 
relates the modes reflected back into section 1 from the junction between sections 1 
and 2 when a mode is incident also from section 1. It is found exactly the same as 
[Rio] except that now the diffraction coefficients are for discontinuities in curvature. 
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For the TE case the matrix elements are given by 

pmn 1. T> T> V 

-**'12 — nI ° 

■ [*W». ».)!«. + (-l) mtn B:(C,»n)l«,] (3.67) 

and for the TM case by 

nmn ^ p D 7 

■**’12 — n “° 

■ [^(ftnAJIe. + (-ir + “CJ(0m,<'n)l9 i ] . (3.68) 

The diffraction coefficients D c th for a discontinuity in curvature are given in Ap- 
pendix C. Notice again that due to the diffraction coefficients being reciprical, [-R12] 
is also symmetric. 

Figure 3.13 shows the relevant geometry for finding [#21]- The matrix [#21] 
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relates the modes reflected back into section 2 from the junction p = 2, q = 1 when 
a mode is also incident from section 2. The derivation is the same as above, except 
that now we are dealing with the modal fields of section 2. It yields the following 
for the TE case, 

R“," = ~\S h ' -R. ” \k m t(p) + h„ a {p)\ - M n }dp (3.69) 

/„ = 0 (3.70) 

M„ = -iMt r {Q t )S{p-k)-iMZ<'(CI«)6(p-a) (3.71) 


h m .(p) = -ijA„Y<,K m (kp) (3.72) 

M*’(Q») = i (3-73) 

V J sm 

«r(Q*M = T-n£Xkb)H"\ki) (3.74) 




1 2 ^ D,(6 am , $an)|Q„ ~opt(s\ \ 

— (y a ) 

J« sm v am 


(3.75) 


^ pt (Qa) 




(3.76) 


~ 2][j2k AmA ’' Y ° 

■ M^^ Kmika)H W {kb)H W(ka) 

Sill v am 

+ Rl(kb)H$(kb)H"Hkb) 
Sin dbm J 


(3.77) 


where the diffraction coefficients are again for a discontinuity in curvature. If either 
of the modes are of the whispering gallery type, one must neglect the second term of 


32 
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Figure 3.14: Double bend geometry in a parallel plate waveguide of infinite extent. 

equation (3.77) which involves the discontinuity at Q\, since there is no whispering 
gallery modal ray incident at Q\,. 

Kt™ for the TM case is found similarly and is given by 

+D k (6.„X«)\QMka)HWXkb)H£\kaj\ (3.T8) 

with the same condition on whispering gallerly modes as before. 

As a numerical example of the reflection from a junction, consider the waveguide 
of Figure 3.14. For the TE case, this guide allows only one propagating mode in 
each section. Therefore, the matrices all have just a single element, so they are 
replaced by scalars. Figure 3.15 plots the magnitude of the component reflection 
coefficients R 12 , R 21 , and iZ 2 4 . Figure 3.16 plots the combined reflection from the 
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Figure 3.15: Reflection coefficients of the three junctions of the double-bend guide 
vs. bend radius. 
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three junctions R ^ SM . It is given by 


r>MSM 

K \2 


R 


MSM 

24 


where 


T » T T>MSM p2 

l 2 i 21^4 r 2 

1 pAf 5 Af p p2 
1 — -ft 24 it2ii 2 

T >2 p p2 
24^21^4 

1 — # 21 ^ 24^4 


Ti2 — \/l — |-Rl2| 2 

r« = y/l - 1-^21 1 2 

*14 = V^ 1 - MW* 
P 2 = P 4 = e“^f. 


(3.79) 

(3.80) 


(3.81) 

(3.82) 

(3.83) 

(3.84) 


The transmission coefficients were found by conserving power flow across a junction. 

As the plots show, these junction reflections are extremely small so that the 
reflection from the termination of an inlet, such as the one in Figure 1.1, would 
completely dominate the total reflection. If the termination is perfectly absorbing, 
the field scattered by just the edges at the open end would dominate total scattering. 
For these reasons, the reflections from the smooth junctions between waveguide 
sections are usually neglected. The amount of energy reflected internally from the 
open end is of the same order as the external scattering by the open end so that 
this scattering mechanism should be included. 


3.2.4 Reflection Matrix Elements for an Impedance Termination 

As an example of a termination reflection matrix, the case of a semi-infinite 
planar dielectic termination is derived in this section. However, the numerical re- 
sults shown later use a perfectly conducting termination because this is a worst 
case where all the energy reaching the termination is reflected back. The reflec- 
tion matrix [i2(] for a planar surface impedance termination is diagonal with the 
elements being the appropriate reflection coefficients for plane waves incident on 
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Figure 3.16: MSM reflection coefficient of the double bend vs. bend radius. 




Figure 3.17: Geometry relevant to finding [Jit]. 


an impedance space. Such an impedance boundary can be created by employing 
a semi-infinite dielectric termination within the waveguide section 5 as shown in 
Figure 1.1. Figure 3.17 shows the geometry relevant to [il t ]. For the TE case the 
diagonal elements are given by 


R? 


cos 0 n ~ \[ibt T ~ sin 2 0 n 


(3.85) 


cos 0 n + ^yjfi r e r - sin 2 0 n 

where fi r and e r are the relative permeability and permittivity, respectively, of the 
termination impedance. For the TM case, the diagonal elements are given by 

cos 0 n - ~ sin 2 0 n 


R? n = 


cos 


0n + ~\fi*r*r - sin 2 0 n 


(3.86) 


It is noted that the surface impedance referred to above is simply that which 
is associated with the above reflection coefficients for the dielectric interface shown 
in Figure 3.17. It is further noted that the expressions for the reflection coefficient 
(and hence the associated impedance) depends on the mode number through 0 n . 
The case of the perfectly conducting planar termination is obtained by letting e r in 
equations (3.85) and (3.86) go to infinity which makes Jl” n equal -1 and 1 for the 
TE and TM cases, respectively. 
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3.2.5 Scattering from the Open End 

Roi describes the scattering from the open end as defined in chapter 2. It is 
found by considering the diffraction of the incident plane wave by the two edges 
of the open end of a semi-infinite parallel plate waveguide. The scattered field, 
including only first order diffractions, can be written as 

U' dge = z[D t>h (ir + 8,n + 8)+ 

+D,, h (ir - 0,* - 9)e-» u *' t ] — (3.87) 

where D, t h is the appropriate soft or hard wedge diffraction coefficient as given in 
Appendix C. Using equations (2.12) and (2.14), the junction reflection coefficient 
can be identified as 

Ro! = D t , h (ir + 6,* + e) + D t>h (‘K-9,*-9)e- i,tkd ''**. (3.88) 
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4. Modal Perturbation Technique 


In this chapter, the hybrid modal analysis method will be extended to inlet 
cavities which have a thin absorber coating on their inner walls. The perturbational 
method which will be used changes only the modal propagation constants /?„ and u n 
and leaves the transverse orthonormal modal field structure unchanged. Therefore, 
the analysis is exactly the same as in the perfectly conducting case of the previous 
chapter, except for small perturbations to the propagation matrices [P,] due to the 
presence of loss on the inlet walls. Also, the diffraction coefficients for the scattering 
from the open end must be adjusted to handle wedges with one perfectly conducting 
face and one impedance face. 

It will be shown for the parallel plate waveguide that the perturbation of the 
propagation constants (3„ can be obtained using only the plane wave reflection co- 
efficients of the impedance surface and the modal ray angle. Extending this idea to 
the annular waveguide makes it easy to obtain the perturbation of the propagation 
constants u n which would otherwise be very complicated both analytically and nu- 
merically. It is then shown that a thin absorber coating can be represented by an 
equivalent surface impedance. 

4.1 Perturbation of the Modal Propagation Constants in a Parallel 
Plate Waveguide Due to Nearly Perfectly Conducting Impedance 
Walls 

Figure 4.1 shows the geometry of a parallel plate waveguide of infinite extent with 
impedance boundary conditions on the walls. Z t is the normalized (to free space 
impedance) surface impedance and Y, — 1 jZ, is the normalized surface admittance. 
Starting with the TE case, the n th orthonormal modal E-field propagating in the 
±*-direction can be written in the form 

e n (x,y) = zP n e n (y)e :f)l3nX (4.1) 
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Figure 4.1: Parallel plate waveguide geometry with impedance walls. 

where P n is a normalization constant. Substituting this into the wave equation 
(A. 4), introducing the separation constant fc„ and solving the resulting differential 
equation for the y variable gives the general solution for the transverse modal field 
as 


e n (y) = e jk ' y + Be~ jk " y (4.2) 

/?« = - kl (4.3) 

B and k n remain to be evaluated using the impedance boundary conditions on the 
walls. The FJ-field must satisfy the equations 

e' n (0)-jkY,e n (0) = 0 (4.4) 

e' n (d) + jkY,e n (d) = 0. (4.5) 


Substituting (4.2) into these boundary requirements gives the solution for B and 
the transcendental equation for finding k n as 


B = 


kn ~ kY, 
k n + kY, 


Y, 2 — 2Y,(-f-) + {~r) 2 = Y, 2 + 2YA + (^) 2 e j2k " d . 


The perturbation requirement is that the surface impedance is nearly perfectly 
conducting, or |F,| 1. Using this approximation, the ( k n /k ) 2 term of (4.7) can 
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(4.8) 


be neglected and (4.2) can be rewritten as 


Y.~ 2(j) = 


Y. + 2I-) 


„}2h n d 


For the perfectly conducting case Y t = oo, the solution for k n is k„ = mr/d. There- 
fore, the perturbed form of k„ will be written as 


— 

" 'fl — 



(4.9) 


where the perturbation k' n is expected to be small for large values of Y,. Substituting 
this form into (4.8) and rearranging gives 


Y, - 2 (kjk) 
Y t +2(k n /k) 


(4.10) 


Because Y, is large, the terms on both sides of (4.10) must be close to unity. There- 
fore, the exponential term can be approximated by the first two terms of its Taylor 
series expansion: 


e i2 < d » 1 + j2k' n d. 


(4.11) 


(4.10) also shows that k' n is small for large Y t , as expected. Using this form in (4.10) 
after some algebra gives 

(j4d/k)k' 2 + (4/fc + j2dY, + j4nir/k)k' n + 4nir/kd = 0. (4*12) 


Again, because Y t is large, all the terms in (4.12) can be neglected except the Y, 
term and the constant term. Solving for k' n gives the perturbation of k n for the TE 
case as 


k' 


j 2 

dY t \kd)' 


(4.13) 


As expected, k' n is small for large values of Y t (it is known from the previous chapter 
that rnr/kd is less than one for propagating modes). 
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The desired result is the perturbation of the propagation constant /?„. This is 
obtained by substituting (4.9) and (4.13) into (4.3) which gives 




(4.14) 


where the (1/ Y,) 2 term is neglected. Because Y, is large, the square root in (4.14) 
can be approximated by the first two terms of its binomial expansion which gives 

j2 (4.15) 


b - Ik* -? 2 4BL- 

P" ~ r { d> 


from which the perturbation to the propagation constant can be identified as 

2 


A = -r 


(S)’ 


Using the modal ray angle expression (A. 17), (3' n can be rewritten as 


(4.16) 


fin = -j 


, 2 sin 2 $ n 


(4.17) 


dY, cos 0 n 

for the TE case. 

The TM case is handled similarly, with the n th orthonormal modal ff -field, 
boundary conditions and transcendental equation for k n represented as follows: 


K(x,y) = zP n h n {y)e* ,<3 "* 

(4.18) 

h n (y) = e ,fc * w + Be~ jknV 

(4.19) 

Pn = sfe - K 

(4.20) 

K(0) - jkZ t h n (0) = 0 

(4.21) 

h' n (d)+jkZ.h n (d) = 0 

(4.22) 

JJ K~ kZ, 

k n + kZ s 

(4.23) 


(4.24) 


The perturbation requirement for this polarization is that the surface impedance 
is nearly perfectly conducting, or Z , is very small. From (4.24) it is seen that Z , 
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should be much less than k n /k. The smallest value of k n is k\ = tt /d (for the 
perfectly conducting case), so the TM perturbation requirement is \Z t \ «C it/kd. 
This is a more difficult requirement to meet because the inlet width d is usually 
fairly large in terms of wavelength. Using this approximation, the Z * term of (4.24) 
can be neglected and (4.24) can be rewritten as 


(j)-2 Z. = 


(j ) + 2Z. 


p j2k n d 


(4.25) 


For the perfectly conducting case Z, = 0, the solution for k n is again k n = mr/d 
and the perturbed form of k n can be written as in (4.9) where the perturbation k' n 
is expected to be small for small values of Z t . It should be mentioned here that 
this perturbation requirement should be heeded much more closely than in the TE 
case because k n can be fairly small. Therefore, Z, should be much smaller than ki. 
Substituting this form into (4.25) and rearranging gives 


(kjjk) 2 Z, _ j2k' n d 
(kn/k) + 2 Z. 


(4.26) 


Because Z t is small, the terms on both sides of (4.26) must be close to unity, so 
(4.11) can again be used to approximate the exponential term. This also shows that 
k' n is small, as expected. Using (4.11) in (4.26) after some algebra gives 


(j2/ k)k'* + ( j2mr/kd + j4Z,)k' n + 4 Z t /d = 0. 


(4.27) 


Because Z, and k' n are small, the k'// and Z,k' n terms of (4.27) can be neglected. 
Doing this and solving for k' n gives 



The perturbation to ft n is found by substituting (4.28) and (4.9) into (4.20) and 
using the binomial approximation to the square root function as before, giving 



(4.29) 
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Using the modal ray expression (A. 17), 0' n can be rewritten as 


& 


. 2 Z._ 

* d cos 0 n 


for the TM case. 


(4.30) 


4.2 Perturbation of the Modal Propagation Constants in a Waveguide 
Using the Plane Wave Reflection Coefficient for an Impedance Sur- 
face 

In this section it will be shown that the perturbation to the propagation con- 
stants can be obtained using only the modal ray angles and plane wave reflection 
coefficients. This will be derived rigorously for the parallel plate waveguide and then 
extended to the more complicated annular waveguide. For this derivation, it is again 
neccessary to assume that the surface impedance is nearly perfectly conducting, i.e., 
IK,! » 1 for the TE case and \Z,\ <C n/kd for the TM case. 


4.2.1 Parallel Plate Waveguide 


Starting with the TE case, the z-propagating term due to the perturbed prop- 
agation constant will be transformed into a term involving only a reflection 

coefficient and the modal ray angle 9 n . The propagating term is first rewritten using 
(4.17) as follows: 




_ 3c *1» 2 9n 
g oi| col 0fi 


g— tm$ n tan^n 


= ( e - * 


3 tan $n 


(4.31) 


Because Y, is large, the term in parenthesis in (4.31) can be approximated with the 
first two terms of its Taylor expansion giving 






| tan 0 n 
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Figure 4.2: Modal ray form of the waveguide modes in a parallel plate waveguide. 


/ - sing n + Y, - ^ sin 2 fl n \ 5 tan *" 

\ sin $ n + Y t ) 

( sin<? n -F t \ SWn 
~ \ sin 0 n + Y») 

= [-i?,(^n)] 5t “' n (4.32) 


where R,(Q n ) is the soft reflection coefficient for a plane wave incident on a surface 
admittance Y, making an angle 9 n with the plane of the surface. The geometrical 
significance of this form is illustrated in Figure 4.2 which shows the crossing plane 
waves of the modal ray form of the waveguide modes of a parallel plate waveguide. 
As the modal field propagates a distance d/ tan 9„ in the ^-direction, the plane waves 
of the modal ray representation undergo one reflection from the waveguide walls. 
This is how the perturbation to the propagation constant of a mode is obtained using 
only a plane wave reflection coefficient and the distance a mode propagates along 
the axis of the guide between consecutive reflections. Using (4.32), this perturbation 
is solved in terms of the reflection coefficient for the TE case as 

Pn = (4.33) 
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(4.34) 


«.(#„) = 


sin 9 n — Y, 


sin 9 n + Y a 

The TM case is done similarly as follows 


g d cos 0 n 


\ tan^n 


_ j»B ilk 

= e d •*** *n 

/ 2Z. 

= I e * in *n J 

M (l 

V sin 9„ J 

' ^K-z.-Sh 

sin 9 n + Z t 

sin9 n - Z a ^' dUn9n 


§ tan B n 


:) 


ft = 3 


\ sin 9 n + z 

= [**(*»)]**“** 
. tan 9 n 


d 


log [R h {9 n )] 


*»(«„) = 


sin 9 n — Z, 


(4.35) 

(4.36) 

(4.37) 


sin 9 n + Z, 

where Rh{9 n ) is the hard reflection coefficient. 

It is mentioned again that the expressions derived in this section for 0' n are valid 
only for the nearly perfectly conducting guide in which |F,| 1 for the TE case 

and \Z t \ < 7r /kd for the TM case. 


4.2.2 Annular Waveguide 

The method derived in the previous section for finding the perturbation to the 
propagation constants using the modal ray angles and reflection coefficients is now 
extended to the annular waveguide. Starting with the whispering gallery type mode, 
Figure 4.3 shows the geometry relevant to finding the perturbation u' n . Considering 
the modal ray form of the modes in this guide, it is seen that the crossing cylindrical 
waves reflect from the outer wall once every angular distance 2 9 an . Therefore, the 
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Figure 4.3: Geometry for finding the perturbation to the propagation constants in 
the annular waveguide, whispering gallery case. 
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perturbed (^-propagation term can be written for the TE case as 


e ^ = [-jR,(0 on )] M «”* . (4.38) 

Solving for u' n gives 

< = lA- log !-«.(#«.)] (4.39) 

*Van 

where the soft reflection coefficient R, is given by (4.34). The TM case is found 
similarly as 

v' n = log («,(«.„)] (4.40) 

“"on 

where the hard reflection coefficient Rh is given by (4.37). 

Figure 4.4 shows the geometry relevant for finding u' n for the non-whispering 
gallery type mode of the annular waveguide. It is seen that a modal ray undergoes 
one reflection each from the inner and outer walls for every angular distance 2(0 an — 
6 bn)- Therefore, the perturbed modal propagation in the ^-direction can be written 
for the TE case as 

e-*** = [R.(0 M )R.(0 h n)}* 9 =£*3 . (4.41) 

Solving for u' n gives for the TE case 

K = 1°8 [*.(«.» (4.42) 

and for the TM case 

K = i , 1 i ( 4 -«) 

AyVan — "bn) 

Once again, it is emphasized that this derivation is valid only for nearly perfectly 
conducting waveguides with |Y,| 1 for the TE case and \Z t \ <C ir/kd for the TM 

case. 
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Figure 4.4: Geometry for finding the perturbation to the propagation constants of 
the annular waveguide, non-whispering gallery case. 
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Figure 4.5: Plane wave reflection from a dieletric layer covering a perfectly conduct- 
ing ground plane. 

4.3 Equivalent Surface Impedance for a Thin Absorber Coating on a 
Perfectly Conducting Ground Plane 

It will be shown in this section that a dielectric layer covering a perfectly con- 
ducting ground plane can be replaced by an equivalent surface impedance using the 
plane wave reflection coefficients for the two types of surfaces [8]. The equivalent 
surface impedance will in general be dependent on the incidence angle. 

Figure 4.5 shows the geometry for this type of configuration with a plane wave 
incident at an angle 0. The soft (TE) reflection coefficient for the coated ground 
plane is given by 

Sin0 — -4-\ArMr — cos 2 0 cot (kty/eTul — cos*0) 

p (0) = >•* - v r ^ r 7 14.44) 

* sin# + — COS* 0 COt {kty/t^ir — COS 2 0) 

where e r and /i r are the relative (complex) permittivity and permeability, respec- 
tively, of the dielectric layer. The corresponding reflection coefficient for a surface 
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impedance is given as in (4.34) as 

R.(0) = 

v ' sin 0 + Y, 


(4.45) 


where Y t = 1 /Z, is the surface admittance. Equating these two types of reflection 
coefficients it is seen that the equivalent surface admittance Y, can be expressed as 
a function of e r , p r , 0 , and kt: 


Y, = - 7 — yc r p f — cos 2 0 cot (kty e r n r — cos 2 6). 

J^T 


(4.46) 


The “thin absorber coating” criterion is fct| v /e r /u r | «C 1. This allows the small 
argument form of the cotangent function to be used which for the TE case gives 


jktfi r 


(4.47) 


which eliminates the angular dependence of the equivalent surface admittance. It 
is seen that for small values of ktfi r , the surface admittance Y, is large, satisfying 
the perturbational requirement of section 4.1. 

For the TM case, the hard reflection coefficients of the dielectric layer covering 
the ground plane and a surface impedance are given by 


r*(«) - 


sin# — ~ cos 2 0 tan (kty/e^fT r — cos 2 0) 

sin0 + — cos 2 0tan (kty/e^JT r — cos 2 6) 


n ( a\ = sin *~ Z ' 
Rh (°> sin 6 + Z.' 


(4.48) 


(4.49) 


respectively, from which it is seen that the equivalent surface impedance can be 
expressed as 


U r n r — cos 2 Otani(ktJe r fi r — cos 2 9). 


(4.50) 


Again using the thin absorber coating criterion, the small argument form of the 
tangent function can be used which gives 


jkt(fi r cos 2 0). 


(4.51) 
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It is seen again that the perturbation requirement of section 4.1 is met if kt^ r is 
small enough to make Z , small. 

For the TM case, Z, is not independent of incidence angle for the thin absorber 
coating case unless e T is large enough to make the second term of (4.51) negligable, 
but still small enough to satisfy the thin absorber coating criterion. If this additional 
approximation is made, then (4.51) is the same as (4.47) because Y, = 1 jZ,. 

4.4 Numerical Results Using the Hybrid 

Asymptotic Modal/Modal Perturbation 
Method 

In this section backscatter results obtained using the hybrid modal technique 
with the perturbed propagation matrices will be presented in the form of backscatter 
vs. aspect angle plots for various inlet configurations. Emphasis will be placed on 
seeing the effects of offsetting the termination of a straight inlet with an S-bend. 
The parameters which will be varied are frequency, polarization, amount of loss 
in the absorber coating and the amount of offset. All of the cases in this section 
have a thin layer of absorber coating on the inner walls of the guide, except for the 
termination which is perfectly conducting. 

Figure 4.6 shows the first two inlet geometries. Notice that the S-bend geometry 
has the same dimensions as the straight inlet with its termination offset one inlet 
width. The inlets will be analyzed at 10 and 35 GHz. At 10 GHz there are 6 
propagating modes in the waveguide sections of these inlets and at 35 GHz there are 
23. Figures 4.7 and 4.8 compare the calculated backscatter vs. aspect angle for the 
lossless straight inlet with measurements taken on a 3-dimensional rectangular inlet 
of the same planar dimensions. The 3-D inlet model was shaped so that scattering 
from the exterior termination was minimized. This scattering mechanism was not 
included in the calculation. The results agree quite well, considering one is in two 
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Figure 4.7: TE backscatter vs. aspect angle for the straight inlet at 10 GHz. — 
Measured, - - Calculated. 



Figure 4.8: TE backscatter vs. angle for the straight inlet at 35 GHz. — Measured, 
- - Calculated. 
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dimensions and the other is in three. 

In the remaining backscatter plots of this section, the perfectly conducting result 
(no absorber coating) will be plotted along with the lossy result and the backscatter 
due to only the open end. The amount of loss due to the absorber coating used is 
approximately 1 dB per reflection at normal incidence. In other words, the power 
density of a plane wave incident on a ground plane coated with this absorber is 
reduced by about 1 dB after reflection. The exact values of the coating parameters 
used here are as follows: At 10 GHz 


t = .001 inches 

e r = 10 -j5 

fir = 5 - jlO 

and at 35 GHz, 

t = .001 inches 
e r = 3 — j3 
fir = 1 )3. 


(4.52) 

(4.53) 

(4.54) 

(4.55) 

(4.56) 

(4.57) 


These values give an equivalent surface admittance of |Y',| ssl7 for the TE case. 
Therefore, the perturbation requirement |F,| >• 1 is fairly well satisfied for these ab- 
sorber parameters. For the TM case, the equivalent surface impedance is \Z,\ «.06. 
ir/kd is .148 at 10 GHz and .042 at 35 GHz. These values do not satisfiy the pertur- 
bation requirement Z, <C ir/kd very well for this case so the TM backscatter results 
for absorber coated ducts in this chapter are not expected to be very accurate. 

Figures 4.9-4.16 show the backscatter vs. aspect angle results for the inlets of 
Figure 4.6 for the two frequencies and polarizations used. The plots show the 
general trend that the effect of the offset is to flatten out the backscatter pattern, 
removing the peak near the axis of the inlet where the termination would be directly 
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Figure 4.9: TE Backscatter vs. aspect angle for the straight inlet at 10 GHz. 
Lossless, 1 dB loss/refl., Rim diffracted. 



Figure 4.10: TE Backscatter vs. aspect angle for the S-bend inlet at 10 GHz. 
Lossless, — 1 dB loss/refl., Rim diffracted. 





Figure 4.11: TM Backscatter vs. aspect angle for the straight inlet at 10 GHz. 
Lossless, — 1 dB loss/refl., - - - Rim diffracted. 



Figure 4.12: TM Backscatter vs. aspect angle for the S-bend inlet at 10 GHz. 
Lossless, — 1 dB loss/refl., - - - Rim diffracted. 
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Figure 4.13: TE Backscatter vs. aspect angle for the straight inlet at 35 GHz. 
Lossless, 1 dB loss/refl., Rim diffracted. 



Figure 4.14: TE Backscatter vs. aspect angle for the S-bend inlet at 35 GHz. 
Lossless, — 1 dB loss/refl., Rim diffracted. 



Figure 4.15: TM Backscatter vs. aspect angle for the straight inlet at 35 GHz. — - 
Lossless, — 1 dB loss/refl., - - - Rim diffracted. 



Figure 4.16: TM Backscatter vs. aspect angle for the S-bend inlet at 35 GHz. — 
Lossless, — 1 dB loss/refl., Rim diffracted. 
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visible for the straight inlet. The presence of the absorber coating further reduces 
the backscatter, however, it does not have much effect in the case of the straight 
inlet near the axis because the perfectly conducting termination is still directly 
visible. 

It is noticed that for the 10 GHz, TM case, the S-bend does not have much 
of an effect on reducing the backscatter near the axis. This is because the TM 
polarization allows the TEM mode (zero mode angle) to propagate in both inlets 
and is the dominant power carrying mode near axial incidence. The S-bend will 
have little effect on this type of mode. Also, the aborber coating will have little 
effect on this mode because the jEMield is normal to the guide walls. However, the 
TEM mode does not exist in more practical three-dimensional inlets, so this is only 
a two-dimensional effect. 

To see if increasing the amount of offset that the S-bend introduces further 
reduces the backscatter pattern, the annular sections of the S-bend inlet of Figure 
4.6 were increased from 42 to 50 degrees in angular extent. The backscatter results 
for this inlet appear in Figures 4.17-4.21. Comparing these results with those of 
Figures 4.10, 4.12, 4.14, and 4.16 it can be seen that increasing the amount of offset 
has little effect on the envelope of the backscatter pattern for the lossless cases 
but reduces further the pattern envelope for the absorber coated cases. This is 
expected because the interior coupled energy has to travel a longer distance inside 
the inlet. For the lossy case this means more attenuation is introduced by the 
absorber coating. However, for the lossless case this only changes the shape of the 
backscatter pattern because the coupled energy is redistributed but none is lost. 
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Figure 4.17: TE Backscatter vs. aspect angle for the extended S-bend inlet at 10 
GHz. — Lossless, — 1 dB loss/refl., - - - Rim diffracted. 



Figure 4.18: TM Backscatter vs. aspect angle for the extended S-bend inlet at 10 
GHz. — Lossless, — 1 dB loss/refl., Rim diffracted. 
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Figure 4.20: TM Backscatter vs. aspect angle for the extended S-bend inlet at 35 
GHz. — Lossless, 1 dB loss/refl., Rixn diffracted. 
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5. Integral Equation Solutions of Internal Wave Reflection and 
Transmission in an S-Shaped Waveguide 

In this chapter, the internal wave reflection and transmission in an S-shaped 
waveguide geometry is solved using an integral equation numerical approach. This 
technique yields a formally rigorous solution, but is not very efficient in terms of 
computer time cost. Nevertheless, it provides a valuable check on the accuracy 
of other approximate and more efficiently obtained solutions such as the one con- 
structed via a hybrid asymptotic modal/modal perturbation approach in Chapters 
3 and 4. 

Referring to Figure 5.1, the S-shaped waveguide consists of two flat and two 
circular sections. The waveguide wall has a normalized surface impedance Z. The 
second flat section may be semi-infinite in length or terminated at a finite distance 
with a normalized surface impedance Z r . The waveguide width is a. The inner 
radius of the circular sections are by and b 3 , and their angular spans are and 
$ 2 * A waveguide mode is incident from the left in the first flat section. Integral 
equations will be set up for solving the aperture fields at the various junctions 
between the different sections. Once the aperture fields at the junctions are known, 
they can be used to calculate the reflected and transmitted modes. 

In the following, presentations will be given for the TM case. But the results 
are also applicable to the TE case, if the impedance Z and Z r are replaced with 
the corresponding admittances. The time convention e? ut is assumed and omitted 
in the formulations. 

5.1 Field Representations for the Flat Sections 

In a parallel-plate waveguide whose wall has a surface impedance Z normalized 
to the free space impedance, the field of a waveguide mode satisfies the homogeneous 
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wave equation 

/ 8 2 d 2 \ 

( dx 2 + dy 2 + fc2 ) Gn ~ ° ~ X ~ a ' ~°° < y < °°» (5.1) 

and the boundary conditions 

dG n 


dx 

dGn 

dx 


— jkZGni &t * — 0, 

= — jkZG n , at * = a, 


(5.2) 


where Jc is the free space wavenumber. The solutions to (5.1) and (5.2) can be 
written in the form 


G n = NnfnWe^V^V, 


where 


/«(*) = 


g-jk/3n* I dn + Z jk0 n x 

Pn ~ Z 


(5.3) 


(5.4) 


i 

I 
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and /?„ are roots of the equation 
2 


(!±!) . . o. 

The waveguide modal field (5.3) is normalized according to 

•>(*)]’*’ = i. 


(5.5) 


such that 


n„ = 


i - s„+zir* 


N)“ 


j2k/3 n a /?„ 

The functions / n (x) associated with different (3 n are orthogonal, i.e. 


(5.6) 


(5.7) 


f fn{x)f m (x)dx = 0, if n ^ m. 

Jo 

A waveguide mode of (5.3) is incident from the left in the first flat section of 
the S-shaped waveguide. The incident wave gives rise to reflected and transmitted 
waves. Therefore, the total field in the first flat section consists of the incident mode 
as well as some reflected modes. 

Let u(xi,yi ) denote the total field (magnetic field in the TM case) where the 
argument (xi,j/i) indicates the first flat section. u(xi,i/i) also satisfies (5.1) and 
(5.2). An expression for u(xi,yi) will be derived with the use of Green’s second 
identity. The desired Green’s function is a solution to the inhomogeneous wave 
equation 

(^2 +~ + fc 2 ) G(*,y ; *',y') = -6(x - x')6(y - y'), 

0 < x < a, 0<y<oo, (5.8) 


satisfying the boundary conditions (5.2) together with dG/dy | tf= o = 0, and the 
radiation condition as |y — y'| — * oo. Specifically, 



- cos( k\J \ - 0^y,)e 


(5.9) 


where y< and y> denote, respectively, the smaller and larger of y' and y. Applying 
Green’s second identity to «(®i,yi) and G{x u yi\ x^y^) in the first flat section yields 


«(*x,!/i) = 2 Nifi(xi) cos(*Vl - 3/1 ) 

VO 

•/.(«» )eW^«», 


_ v jg. 

n jka 


dy'i 


y[=o 


dx j 


(5.10) 


where the subscript i is the index number of the incident mode. 

The same procedure can be readily applied to obtain an expression for the total 
field u(x 2 ,y 2 ) in the second flat section, 


«(**»!&) = 


N* 1 — 


jka 1 + Rn 
i \ 9u{x' 2 ,y' 2 ) 


/* /»(«',) 

v0 


9y 2 


dx' 


y'i=o 


(5.11) 


where i£„ = 0 if the second flat section is semi-infinite, and 


= Vv(EjL -W^, 

z. + s/i- K 


if the second flat section is terminated with a surface impedance Z r (normalized to 
the free space impedance) at a distance r from the junction between the second flat 
and curved sections. 
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and the boundary conditions 

3G 


dp 

6G 

dp 

dG 

d<t> 


= jkZG, at p — b, 

= —jkZG, at p = b + a, 


= 0, at <f> — 0, and <f> — $. 


The solution for G(p , <j>; p', <f>') is 

G(p,<j>\p\<j>') = X) *^oP>)cos(^)cos(t^), 

where 

Pn „ 

v — P = an integer, 

$ 

, x _ [g^>(feg<) - BJOJep^)\ [CyJ v ( kp>) ~ ^i\r w (fe/g>)] 
9v W<' p> > eov(A„D v - B V C„) 

J v , N w = Bessel and Neuman functions, 

K, N' v = derivatives of J„ and N„ w.r.t. the argument, 

P<iP> — the smaller and larger of p and p\ 

A v = J' v {kb + ka) + jZJ v (kb + ka) 

B v = J' v (kb)-jZJ v {kb) 

G v — Nl(kb + ka) + jZJ v {kb + ka) 


( 5 . 13 ) 


( 5 . 14 ) 


D v = N' v (kb) - jZN„(kb) 

= 0 

0 ’ 




* / J ~ * 

_ j 2, if v = 

— \ 1, if v > 


Now, applying Green’s second identity to G(p,4>-,p' ,<(>') and u(p,<f>), which sat- 
isfies the homogeneous wave equation 

d 2 


1 d_ ( d\ 1 
pdp v dp) + p‘ 


+ k 2 


2 d<j > 2 

0<<f><Q,b<p<b + a, 


u(p,<f>) = 0, 


( 5 . 15 ) 
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and the boundary conditions 


du 

dp 

du 

dp 


jkZu, at p = b, 

— jkZu , at p = b + a, 


yields an expression for u(p,<j>) y 


7T 

u(p,<j>) = X, $ cos(i/<^) j b g u (p < , 


V>>) 


(-ir 


p'd<f>' 


<=* 


p'd# 


dp'. 


I^'=0j 


(5.16) 


(5.17) 


One may attach a subscript 1 or 2 to the various variables in (5.17) to signify either 
one of the two circular sections. 


5.3 Integral Equations for the Unknown Aperture Fields 


One may set up integral equations for solving the unknown aperture field quanti- 
ties du/dy (or du/d<t>) at the junctions from (5.10), (5.11), and (5.17), by imposing 
the following boundary conditions 


0«(*i>yi) 


dy x 

du{x 2 ,y 2 ) 


du{pi,<t> x ) 


d y 2 

du(p u <f>i) 


Pid<j>i 


VI =0 


V2=0 


4>i=0 


P\d<i>\ 

du{p 2 ,<}> 2 ) 


= Mi(x x ) = M x (px), 


<h=*i 


Ptdfa 

< MptM 

Pid<j > 2 


= —M 2 (x 2 ) = —M 2 (p 2 )t 


<h=*2 


= Ms(pi) = M 3 {p 3 ), 


(5.18) 


<h=0 


and 


«(*i,yi)lvi=o = «(/>i,<Ai)k=*i» 
«(*2,y2)| W =0 = u(p 2 ,<f> 2 )\<h=* 7 i 

tt(/»li^l)U =0 = tt(/> 2 i^a)l*i= 0 - 


(5.19) 
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The integral equations are 

+ E^J^ a 9Ap<,P>)[M 1 (p[)-(-l) p M 3 (p[)]dp'i = 2 Niffa), (5.20) 

£ L 9v{p< ,?>)[(- 1 ) P Mx (p[ ) - M 3 (p [ )] dp[ 

+ }AP<,P>)[(-l) r MM)- mp$dpl = 0, (5.21) 


? jSfr x Mm,) Io 

+ £ £/“%..(/><, I>>) [M,W) - (-lJ'XW)] dp’, = 0. (5.22) 

Note that for simplicity, the subscripts 1 and 2, which distinguish terms in the two 
circular sections, are omitted from some of the variables under the summation 
in (5.20)-(5.22). 

The integral equations (5.20) - (5.22) will be solved by the moment method. 
Using the substitutions 

MiM = f)PJk/.(*i) 

« 

mm) = £cw.(f>i-») 

n 

MM) = JW.(«,) 

n 

MM) = Y.vmm-1) 

n 

MM) = £ W n kfM ~ <■) 

n 

JV»0>j) = f)WJkA(« + »-/>,), (5.23) 

n 

where U n , U n and Wi, are unknown constants, multiplying (5.20)-(5.22) by kf m (x), 
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and integrating the resulting equations from * = 0 to * = a, one obtains 

*° U m + £; - w „ r ° mn ) = , 2t< * 






E - W'nV.U) + E < 

n Pm A n 


l-Rm 


ka 


i + a . j ( i -«,)«2 


v „ + E «'**•*) 


m = 1,2, ... ,iV , 


(5.24) 

0, 

(5.25) 

0, 

(5.26) 


where 

Iqmn = J o kdxf m (x)J o kdx'f n (x , )'^2^-9u(b + z < ,b + x > ), 

Yqmn = J kdxf m (x)J^ kdx’ f^x’)^ ^-(-l) p g„(b + x<,b + x^, 

q = l,2. (5.27) 


In sum, there are 37V linear equations which can be solved for the unknown U n , V n 
and W n . 

Using the substitutions given in (5.23), the reflected waves in the first flat section 
can now be written as 

«’(*■, yi) = 

- E ; D \ „ *■./»(*■ (5.28) 

N umerical examples of the reflection coefficients associated with (5.28) are given 
next. 


5.4 Numerical Results 

The integral equation technique presented in this chapter is not efficient in terms 
of the computer time cost. The primary reason is due to the fact that the series 




in the field representation for the circular section given in (5.17) converges very 
slowly 1 . For instance, in the following examples, 12 terms are retained in the series 
representation for the fields in the flat sections given in (5.10) and (5.11), in contrast 
to about 200 terms retained in the series of (5.17). Nevertheless, as pointed out 
earlier, this formally exact solution is presented mostly to check the accuracy of the 
perturbed hybrid approach of Chapter 4. 

Table 1 lists the reflection coefficients Rn m for the waveguide geometry speci- 
fications: a = 3.386667A, = 5, = 4.766733A, = $ 2 = 45°, Y = 5.046133A, 

Z = Z r = 0, TM case. MM denotes the moment method solution of this chapter 
and HM denotes the hybrid solution of chapters 2, 3 and 4. 

1 The computer program for computing Bessel and Neuman functions of fractional orders was 
supplied by Prof. J.H. Richmond. 
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TABLE 1 


nm 

Rnm (MM) 

Rrtm (HM) 

i 

1 

0.644 Z45.8® 

0.649 Z42.9® 

i 

2 

0.253 Z68.9® 

0.246 Z65.4® 

i 

3 

0.262 Z-144.4® 

0.257 Z-148.0 1 

i 

4 

0.576 Z-20.4® 

0.578 Z-23.6® 

i 

5 

0.095 /.-59.6 0 

0.099 Z-59.6® 

i 

6 

0.280 Z-0.6® 

0.283 Z-3.2® 

i 

7 

0.182 Z-16.0® 

0.172 Z-16.7® 

2 

2 

0.880 Z-173.6 0 

0.885 Z-177.1' 

2 

3 

0.158 Z-l 16.0° 

0.141 Z114.9® 

2 

4 

0.218 Z-117.2 0 

0.217 Z-120.6* 

2 

5 

0.054 Z146.3 0 

0.056 Z141.6® 

2 

6 

0.209 Z-100.0® 

0.204 Z-103.9' 

2 

7 

0.206 Z-123.9 0 

0.217 Z-127.0' 

3 

3 

0.658 Z175.5 0 

0.660 Z172.2® 

3 

4 

0.187 Z179.1® 

0.187 Z178.3® 

3 

5 

0.571 Z- 3.7° 

0.564 Z-6.9® 

3 

6 

0.342 Z37.3° 

0.347 Z35.4® 

3 

7 

0.088 Z-96.3 0 

0.083 Z-98.7® 

4 

4 

0.725 Z92.7® 

0.723 Z90.0® 

4 

5 

0.151 Z 179.7° 

0.155 Z 177.0® 

4 

6 

0.150 Z17.2° 

0.158 Z17.2® 

4 

7 

0.119 Z-110.1® 

0.117 Z-115.0' 

5 

5 

0.383 Z37.6® 

0.392 Z34.3® 

5 

6 

0.691 Z-10.6® 

0.684 Z-12.6® 

5 

7 

0.165 Z-56.5® 

0.153 Z-59.2® 

6 

6 

0.472 Z166.3® 

0.486 Z164.7® 

6 

7 

0.197 Z2.1® 

0.182 Z-0.9® 

7 

7 

0.915 Z107.1® 

0.914 Z 104.3° 


The reflection coefficients for the same waveguide geometry but with the TE 
polarization are listed in Table 2. 
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TABLE 2 


nm 

7?nm 

, (MM) 

Rnm (HM) 

i 

1 

0.722 

Z51.8 8 

0.716 Z50.6 8 

i 

2 

0.296 

Z94.1 8 

0.288 Z94.4 8 

i 

3 

0.409 

Z96.7 8 

0.505 Z92.6 8 

i 

4 

0.237 

Z-135.3 0 

0.248 Z-138.0 8 

i 

5 

0.266 

Z96.9 0 

0.265 Z91.7 8 

i 

6 

0.127 

Zl 13.8° 

0.130 Z109.9 8 

2 

2 

0.552 

Z-109.4 0 

0.543Z115.0 8 

2 

3 

0.339 

Z-24° 

0.355 Z-29.9 8 

2 

4 

0.591 

Z-123.1 0 

0.586 Z-126.3 8 

2 

5 

0.337 

Z-21 .5° 

0.349 Z-25.1* 

2 

6 

0.174 

Z-23.3 0 

0.171 Z-27.9 8 

3 

3 

0.515 

Z-47.3 0 

0.489 Z-52.6 8 

3 

4 

0.435 

Z95.9 8 

0.426 Z92.4 8 

3 

5 

0.368 

h-* 

o* 

00 

o 

0.372Z155.5 8 

3 

6 

0.219 

Z64.9 8 

0.222 Z61.8 8 

4 

4 

0.477 

Z119.1 8 

0.501 Z118.3 8 

4 

5 

0.384 

Z93.1 8 

0.359 Z87.2 8 

4 

6 

0.177 

Z91.4 8 

0.179 Z90.6 8 

5 

5 

0.723 

Z-36.3 8 

0.727Z-40.6 8 

5 

6 

0.104 

Z116.5 8 

0.101Z113.9 8 

6 

6 

0.929 

Z-88.1 8 

0.927Z-90.8 8 


Two examples of the reflection coefficients for the waveguide geometry with 
impedance wall are given in Tables 3 and 4. The geometry specifications are: a = 
3.386667A, bi = b 2 = 4.766733A, = * 2 = 42°, Y = 5.842A, Z r = 0, Z = 

0.0532 + j‘0.0266, TM case in Table 3 and TE case in Table 4. 
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TABLE 3 


nm 

Rnm (MM) 

Rnm (HM) 

i 

1 

0.324 Z49.0° 

0.319 Z-17.4 0 

i 

2 

0.113 /-35.8° 

0.088 Z-58.2 0 

i 

3 

0.064 Z25.4° 

0.044 Z-72.9 0 

i 

4 

0.057 Z-143.9 0 

0.058 Z-119.7 0 

i 

5 

0.130 Z118.6 0 

0.115 Z110.6° 

i 

6 

0.020 Z66.5 0 

0.029 Z64.9° 

i 

7 

0.007 Z-79.9 0 

0.006Z-62.4 0 

2 

2 

0.012 Z 169.7° 

0.049/150.5° 

2 

3 

0.303 Z44.7 0 

0.284 Z43.8° 

2 

4 

0.042 Z-163.9 0 

0.041Z119.5 0 

2 

5 

0.178 Z88.2° 

0.191 Z79.3° 

2 

6 

0.067 Z98.9 0 

0.076 /89.1° 

2 

7 

0.007 Z102.8 0 

0.013/85.9° 

3 

3 

0.249 /67.4° 

0.267 /55.3° 

3 

4 

0.108 Z-141.1 0 

0.49/453.5° 

3 

5 

0.194 Z-140.3° 

0.179/440.6° 

3 

6 

0.019 / 155.3° 

0.041 /-166.6 0 

3 

7 

0.028 Z-0.8° 

0.026 Z- 2.2° 

4 

4 

0.385 Z7.8° 

0.406 Z6.0° 

4 

5 

0.104 Z-172.5 0 

0.093/475.8° 

4 

6 

0.008 Z130.7° 

0.049 Z462.0° 

4 

7 

0.004 Z69.3° 

0.013 Z47.9° 

5 

5 

0.130 Z64.2° 

0.127 /79.6° 

5 

6 

0.190 Z175.2 0 

0.269/172.6° 

5 

7 

0.055 Z89.4° 

0.061 Z84.4° 

6 

6 

0.221 Z57.6° 

0.329 Z61.9 0 

6 

7 

0.012 Z-71.7 0 

0.015/-88.5 0 

7 

7 

0.144 Z13.1° 

0.137/11.6° 
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TABLE 4 


nm 

Rnm 

(MM) 

Rnm (HM) 

i 

1 

0.450 

Z39.0® 

0.479 Z42.7° 

i 

2 

0.120 

Z-74.7® 

0.149 Z-68.5® 

i 

3 

0.368 

Z57.4® 

0.364 Z54.6 0 

i 

4 

0.120 

Z 154.2° 

0.105 Z 152.5° 

i 

5 

0.065 

Z-12.3 0 

0.074 Z-6.0° 

i 

6 

0.028 

Z-153.7® 

0.016Z162.8 0 

2 

2 

0.451 

Z-150.4 0 

0.476Z-146.7 0 

2 

3 

0.288 

Z148.2 0 

0.264 Z152.7 0 

2 

4 

0.370 

Z160.7° 

0.386 Z155.7® 

2 

5 

0.197 

Z-84.0 0 

0.195 Z-86.2® 

2 

6 

0.089 

Z-174.3® 

0.080Z-172.5® 

3 

3 

0.612 

Z-127.5 0 

0.584Z-130.4® 

3 

4 

0.234 

Zll.l 0 

0.243 Z18.5® 

3 

5 

0.052 

Z-18.4 0 

0.057Z-36.2® 

3 

6 

0.036 

Z-138.2 0 

0.030Z-140.2® 

4 

4 

0.569 

Z- 72.6° 

0.583 Z-69.1® 

4 

5 

0.153 

Z-42.9 0 

0.184 Z-52.4® 

4 

6 

0.078 

Z-85.8 0 

0.067 Z-83.9® 

5 

5 

0.705 

Z171.2 0 

0.707Z171.3® 

5 

6 

0.088 

Z38.6 0 

0.075 Z37.4® 

6 

6 

0.249 

Z-165.6 0 

0.238Z-159.2® 


Note that for the TE case, the appropriate corresponding admittance must be 
substituted for Z and Z T . As the above comparisons show, the more efficient hybrid 
asymptotic modal/modal perturbation method closely approximates the formally 
exact moment method solution. 
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0. Alternative Solution for the Scattered Field Using the Geometrical 

Optics Ray Approach 


In this chapter the electromagnetic backscatter from the inlet of Figure 1.1 
is calculated using the geometrical optics ray technique. The geometrical optics 
(GO) field incident on the inlet opening is partly scattered in the exterior region 
by the edges at the open end as U edae and the rest is coupled into the inlet. Only 
the incident rays which directly enter into the inlet are included; the remaining 
contribution to the fields coupled into the inlet via diffraction of the incident field 
by the edges at the open end is ignored. The incident rays which enter into the 
inlet undergo multiple reflections from the interior walls of the inlet. These GO rays 
which are multiply reflected from the interior walls of the inlet propagate toward the 
termination at the rear of the inlet and then bounce back via additional reflections 
at the walls to provide a GO based field representation in the aperture at the open 
end. Further wave interactions between the open front end and the rear termination 
are ignored because these are weaker. The field U eav which is essentially produced 
by the radiation into the exterior space from the above mentioned GO based field 
in the aperture at the open end, is found as usual via the aperture integration (AI) 
technique [11,12]. This aperture integration is done exactly the same way as in 
chapter 3 using equation (3.47). The difference is that now U ap is based on the GO 
field in the aperture at the open end rather than the modal field of equation (3.48). 
The advantage of this method is that although not as accurate as the other methods 
discussed earlier, it can handle inlets of arbitrary shape and absorber coating, rather 
than ones that are made up of uniform waveguide sections and have thin absorber 
coatings. Also, as mentioned breifly in chaper 2, the hybrid asymptotic modal 
method becomes increasingly cumbersome as the guide width becomes large. The 
GO/AI method actually improves in accuracy for those cases. On the other hand, 
the GO/AI method can also become inefficient and cumbersome for inlets because 


76 



of the large number of rays to be tracked. The main reason for any inaccuracies in 
this GO/AI approach is because it does not include any effects of diffraction that 
enter within the duct and other higher order effects. The latter include the junction 
and aperture reflections considered in chapter 2. The GO/AI method of analysis is 
presented here for comparison purposes and to supply an alternative procedure for 
analyzing the scattering from inlet cavities, especially at higher frequencies where 
the modal method may become too inefficient and for more lossy absorber coatings 
for which the modal perturbation technique is not valid. Again, U stands for the 
z- component of the 25-field for the TE case and the z-component of the 22-field for 
the TM case. 

6.1 Tracking the GO Ray Field 

The GO field is found by breaking up the incident plane wave which enters the 
inlet into N collimated beams of equal width, as shown in Figure 6.1. These smaller 
columnated beams or ray tubes are traced through the inlet via the rear termination 
until they are about to exit through the aperture, as shown in Figure 6.2. Actually, 
only the two rays which form the boundaries of the smaller ray tubes are traced. 
The contribution of each ray tube to the aperture integral is found by breaking up 
the AI equation of Appendix E into a summation of sub-apertures formed by each 
of the exiting ray tubes in the aperture. 

U cav = J 3 -cos9—-Y, U?(y)e- jkv ' m0 dy (6.1) 

V 2lr VP n=l 

where each integral in the summation represents the contribution to the field scat- 
tered by the interior cavity effects arising from each ray tube. It is easy to show 
that the two terms of (E.6) and (E.12) are the same for this formulation which is 
why (6.1) has only one term. Figure 6.3 shows the relevant geometry pertaining to 
the limits y n and y n+ 1 which occur in the integrals of this summation. The y n and 
y n+ 1 shown are found from ray-tracing. Notice that because the ray tube undergoes 
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Figure 6.1: Incident plane wave divided into N columnated beams of width A / . 



Figure 6.2: Ray tube traced through an inlet. 

! 

I 
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reflections from curved inlet walls, it is no longer columnated and has some spread- 
ing (either divergent or convergent). The field in the ray tube will appear to be due 
to a radiating line source located at some caustic, as shown, with a phase factor 
representing the propagation path delay of the ray tube up to that point. The mag- 
nitude of the field in each subaperture (y n < y < y„+i) is obtained by conserving 
power in the ray tube which terminates in that subaperture. It is assumed that 
each ray tube is small enough so that the magnitude of the field in the sub-aperture 
is approximately constant within that subaperture. To conserve power in each ray 
tube, 

( 6 - 2 ) 


where A' is the width of the ray tube when it enters the inlet, and it is given by 

A' = 4 cost?; (6.3) 

N 

whereas, A n is the width when it exits, transverse to the direction of propagation. 
U* p (y) is the scalar part of U“ p (y). The (scalar) field in the aperture is then given 
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via GO by 


V?(y) 



P n jCn e -jk[t(y)+Rn] 


(6.4) 


where R n is the path length of the n th ray tube inside the inlet up to the last caustic, 
r n is the product of all the reflection coefficients which describe the reflections 
encountered each time the ray tube reflected from an inner wall of the inlet, and c n 
is the total number of caustics the n th ray tube passed through. The rays are traced 
according to Snell’s law, which requires that ray paths before and after reflections 
are straight lines, and at reflection points the angle of reflection equals the angle of 
incidence. Furthermore, the reflected ray, the surface normal and the incident ray 
at the point of reflection all lie in the same plane which for the two-dimensional case 
is always in the plane of the paper. The reflection coefficient for a ray reflecting off 
a semi-infinite dielectric termination inside the inlet cavity is given for the TE case 
by 


T t 


cos<t> - ^ - sin 2 <f> 

c°s<t> + u r e r - sin 2 <f> 


and for the TM case by 


(6.5) 


cos6 — —Ja r t T — sin 2 6 

Tt ? V ■ (6 . 6 ) 

cos<f> + ^ \J n r e r — sin <f> 

where jU r and e r are the relative permeability and permittivity of the impedance, 
respectively, and <f> is the angle the n th ray tube makes with the surface normal 
when it strikes the dielectric termination. The reflection coefficient for a dielectric 
coating on a ground plane is given by [8] 


sin# — — cos 2 # cot (kiy/e^fTr — COS 2 #) 

sin # + yj ~€rjxf — cos 2 # cot (kt\Ze^/M r — cos 2 #) 


(6.7) 
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for the TE case and by 

^ sin# — ^-v^r^r — cos 2 # t<a\{kt\/e r [i T — cos 2 #) 

sin# + ^■v'e r /t P — cos 2 # ta.n(kt^/e^T r — cos 2 #) 

for the TM case where e r and p r are the relative permittivity and permeability, 
respectively, of the coating, t is the thickness of the coating, k is the free space wave 
number and # is the angle the propagation vector of the incident plane wave makes 
with the surface. 

To find c n , the total number of caustics a ray passes through, it is necessary to 
keep a track of the relative caustic position as the ray is traced. Each reflection 
from a curved surface changes this position. Figure 6.4 shows a typical reflection 
from such an interior boundary or surface of the inlet and how a new caustic is 
formed. The new caustic distance p e after that reflection in Figure 6.4 is given in 
terms of the previous caustic distance p' e before the reflection in that figure by [5] 

_ p' c R sin <f> 

R sin <f> + 2 p' c 

Notice that p e as well as p' c can be negative. R which is the radius of curvature 
of the surface at the point of reflection is positive for a convex boundary but is 
negative for a concave boundary. If p e < 0 and \p c \ is less than the distance from 
the reflection point to the reciever point p, then the caustic lies inside of the inlet 
and c n is incremented. This must be checked after each reflection while tracing the 
ray through the inlet. The caustic of the incident ray tube entering the inlet is at 
infinity in the direction of the source. 

The integral in equation (6.1) can now be written as: 




f Vn+1 U“ p {y)e- jkvs ' m0 dy 

J Vn 



p ^jCn g—jk(R n -l n ) . 





( 6 . 10 ) 
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where 


s(y) = \f% + (y - y n ) 2 + 2(y - y„)s n sin0 n . 

For Ay n <C s n , s(y) can be approximated as a linear function of y 
*{y) « «n + (y - y n ) sin 0 n 


( 6 . 11 ) 


( 6 . 12 ) 


so that the integral of equation (6.1) can be evaluated in closed form to give the 
final result as: 


f Vn+l U° p (y)e~i kv * mt dy = p n j c " C -i fe [- B «+|(v»+vn+>) eintf + j Ay B tin*„] # 

'vn n V A n 

•At/ n sinc ^A:Ay„(sin^„ + sin^)j . (6.13) 

The total scattered field is given as the sum of the cavity scattering and the field 
scattered by the edges at the open end of the inlet: 


•jjBcat _ ^jcav _j_ jjcdge 


U edae = zRo^* 




(6.14) 

(6.15) 


where Rqi is given by equation (3.88). 


6.2 Numerical Examples Based on the GO/AI Method 

As it turns out, the GO/AI method gives increasingly less accurate results as 
the inlets become longer. This is due to the fact that the method does not include 
the diffracted fields inside the inlet which generally become dominant as the wave 
propagation distance inside the inlet increases. The reflected beam that is initially 
collimated usually becomes diffused at larger distances inside the inlet because of 
the presence of fields diffracted by edges at the open end which enter into the 
inlet. Also, the method becomes less accurate at lower frequencies because fewer 
inlet waveguide modes propagate. For these reasons, the sample inlets used in this 
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Figure 6.6: Offset inlet geometry. 

section are shallower than those of the previous chapters. Their geometries are 
shown in Figures 6.5 and 6.6; it is seen from those figures that again one inlet is 
straight and the other is offset or S-shaped in which the offset inlet consists of two 
curved sections. The parameters of the absorber coating are the same as in (4.52)- 
(4.57) with an additional case where t = .005”. This lossier coating causes about 
5 dB of loss per reflection at normal incidence compared with 1 dB for the thinner 
coating ( t = .001”). 

Figures 6.7-6.22 show the backscatter pattern vs. aspect angle for the two inlets 
of Figures 6.5 and 6.6 at 10 and 35 GHz for the two polarizations. The top figure 
on each page is the hybrid solution for the lossless and 1 dB loss/reflection cases, 
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included for comparison. The aspect angle is limited to within 20 degrees 

of the axis and the data points are 2 degrees apart in these calculations because 
the ray-tracing method takes a large amount of computing time. In general, the 
more complex the inlet shape, the more ray tubes are needed to converge to the 
GO based result. Therefore, the offset inlet of Figure 6.6 requires more ray tubes 
(and hence more computing time) than the simple parallel plate inlet of Figure 6.5. 

The plots illustrate that the GO/AI method can roughly predict the average 
of the backscatter pattern, but details are lacking. This is expected due to the 
approximate nature of a technique which does not include effects of diffraction. 
The only advantage of this method, other than its analytic simplicity, is that it is 
not limited to inlets made up of uniform waveguide sections which a modal method 
would require. It is noticed from the backscatter patterns that for this shallower 
inlet, the termination is not as well hidden for near-axial incidence and that the 
presence of the absorber coating has a substantial effect in consistently reducing 
the overall backscatter pattern. 


85 




Figure 6.7: TE Backscatter pattern for the straight inlet at 10 GHz, modal solution. 
— Lossless, — 1 dB loss/refl., - - - leading edges only. 



Figure 6.8: TE Backscatter pattern for the straight inlet at 10 GHz, GO/AI so- 
lution. — Lossless, 1 dB loss/refl., - - - 5 dB loss/refl., • • • leading edges 

only. 
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Figure 6.9: TE Backscatter pattern for the S-bend inlet at 10 GHz, modal solution. 
— Lossless, — 1 dB loss/refl., - - - leading edges only. 



Figure 6.10: TE Backscatter pattern for the S-bend inlet at 10 GHz, GO/AI so- 
lution. — Lossless, 1 dB loss/refl., 5 dB loss/refl., • • • leading edges 

only. 
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Figure 6.11: TM Backscatter pattern for the straight inlet at 10 GHz, modal solu- 
tion. — Lossless, — 1 dB loss/refl., leading edges only. 
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Figure 6.12: TM Backscatter pattern for the straight inlet at 10 GHz, GO/AI 
solution. — Lossless, — 1 dB loss/refl., - - - 5 dB loss/refl., ••• leading edges only. 
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4b. -it. -n. -*. o. *. i«. ii. n. 

THETH (DEG) 

Figure 6.13: TM Backscatter pattern for the S-bend inlet at 10 GHz, modal solution. 
— Lossless, — 1 dB loss/refl., leading edges only. 



Figure 6.14: TM Backscatter pattern for the S-bend inlet at 10 GHz, GO/AI so- 
lution. — Lossless, 1 dB loss/refl., 5 dB loss/refl., • • • leading edges 

only. 
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Figure 6.15: TE Backscatter pattern for the straight inlet at 35 GHz, modal solu- 
tion. — Lossless, — 1 dB loss/refl., leading edges only. 
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Figure 6.16: TE Backscatter pattern for the straight inlet at 35 GHz, GO/AI 
solution. — Lossless, — 1 dB loss/refl., - - - 5 dB loss/refl., • • • leading edges only. 
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Figure 6.17: TE Backscatter pattern for the S-bend inlet at 35 GHz, modal solution. 
— Lossless, — 1 dB loss/refl., - - - leading edges only. 



Figure 6.18: TE Backscatter pattern for the S-bend inlet at 35 GHz, GO/AI so- 
lution. — Lossless, 1 dB loss/refl., - - - 5 dB loss/refl., • • • leading edges 

only. 
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Figure 6.19: TM Backscatter pattern for the straight inlet at 35 GHz, modal solu- 
tion. — Lossless, — 1 dB loss/refl., - - - leading edges only. 



Figure 6.20: TM Backscatter pattern for the straight inlet at 35 GHz, GO/AI 
solution. — Lossless, — 1 dB loss/refl., 5 dB loss/refl., ••• leading edges only. 
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Figure 6.21: TM Backscatter pattern for the S-bend inlet at 35 GHz, modal solution. 
— Lossless, — 1 dB loss/refl., - - - leading edges only. 



Figure 6.22: TM Backscatter pattern for the S-bend inlet at 35 GHz, GO/AI so- 
lution. — Lossless, 1 dB loss/refl., 5 dB loss/refl., • • • leading edges 

only. 
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7. Conclusion 


It has been shown that the effect of employing an S-shaped inlet for making 
the termination not directly visible alters the backscattered pattern so that it does 
not peak near the forward axis of the inlet. Also, the overall pattern envelope 
flattened out and was reduced by a few dB, especially for longer inlets. Adding 
the absorber coating to the straight inlet did not effectively remove the peak in the 
backscatter pattern near the inlet axis. Combining the absorber coating with the 
S-bend offset effectively and consistently resulted in flattening out and reducing the 
overall backscatter pattern, even near axial incidence. 

The scattering from the interior inlet cavity termination is the dominant con- 
tributor to the total scattering by the inlet cavity except when the termination is 
highly absorbing in which case the scattering from the edges at the opening in front 
of the inlet dominates. It was also shown that the reflections from the junctions of 
waveguide sections from which the S-shaped inlet is “built up” are very small and 
can usually be neglected. The interior reflection from the open end was seen to be 
on the order of the scattering by the edges so this contribution to the backscatter 
should be included. 

The field expansion/matching technique of chapter 5 gives analytically accurate 
results and can handle any absorber coating which can be characterized by an equiv- 
alent surface impedance. However, the method is largely numerical and therefore 
gives little insight into the scattering mechanisms of the interior region of the inlet 
as well as requiring a large amount of computing time. Its main purpose here was 
as a means of checking the more approximate and efficient hybrid technique. 

The GO/AI method of finding the backscattered field was not as accurate as the 
modal method, especially for smaller guide widths and longer inlets. However, its 
usefulness becomes apparent when the guide width becomes large so that there are 
many propagating modes. At higher frequencies, the modal method becomes more 
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cumbersome but it remains quite accurate. The GO/AI method tends to become 
more accurate at higher frequencies. It is reasonable to suggest that the modal 
method be used for guides of width less than about 15 wavelengths (30 propagating 
modes) and the GO/AI method for larger guide widths. 
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A. Orthonormal Modes of the Parallel Plate Waveguide 


This appendix derives the expressions for the orthonormal modal fields of the 
uniform parallel plate waveguide. The modal ray form of these modes is also derived. 

“Orthonormal” means that each mode carries power independently of all the 
other modes and they are normalized to carry unit power. Mathematically, this is 
written 

I e m (y ) x h n (y) • dS = S mn (A.l) 

Js 


where S mn 

6 = 


is the Kroeneker delta function given by 

{ 1 if m = n 
0 if m / n 


(A.2) 


and S represents any complete cross-section of the guide. 

Figure A.l shows a parallel plate waveguide of infinite extent. Starting with the 
TE case, the electric field is in the z-direction and is a function of x and y 


E(x,y) = zE,(z,y). 


(A.3) 


The FJ-field is a solution to the wave equation 
(V 2 + k 2 )E(x,y) = 0 


(A.4) 
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Figure A.l: Parallel plate waveguide geometry. 
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which in scalar form is 


(£5 + h + * 2 ) E,( * ,v) = 0 (A5) 

where k is the free space wave number, 2ir/\. The boundary condition is that E z 
equals zero on the walls of the guide. Using separation of variables and enforcing 
the boundary conditions the solution is found to be an infinite summation of modes 

E(*,y) = EUteMe-’^ + A-U (A.6) 

n=l 

where 


e n (y ) = £P„ sin(— t/) (A. 7) 

A. = yjk 2 -^) 1 - (A.S) 

P n is a constant which normalizes the power of the n th mode and A* are arbitrary 
constants. The infinite sum is usually truncated to include only propagating modes. 
For n greater than some number N , the modes are evanescent in nature and die 
out exponentially along the axis of the guide. This cutoff number is one less than 
the the value of n which makes /?„ pure imaginary, or 
kd 

N = int(— ) (A.9) 

7T 

(the integer portion of kd/ir). 

The H -field of the guide is derived directly from the .E-field using one of Maxwell’s 
equations 


V x E(x,y) 


A 

Y 0 


which gives 


H(x,y) 


n 

jk 



(A.10) 
(A. 11) 
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(A.12) 


= £ [- 4 A(!') e + 4.M») e “’"'] 


n = 1 


&i(») = an 


'.Pn ■ . n * , . »* - 

t »T”( 7 ! ' ) - I iW C05( T ! ' ) 


(A. 13) 


P n is found by normalizing to unity the power carried by the n tH mode. This is 
done by integrating the Pointing vector over a cross-section of the guide and setting 
this equal to unity, 

r d 

/o 

Substituting and integrating gives 


/ e„(y) x h±(y) • (±*)dy = 1. 
Jo 


(A.14) 


Pn = 


!2kZo 

~W' 


(A.15) 


The ray-optical form of these modes is found by writing the sine function of 
equation (A.6) in its Euler (exponential) form and combining the exponentials, 


e n (y)e* jl3nX = z~ le^rv^nm) _ e j(-^vT0nx)] 

2 j L J 


By making the substitution 

nit 


sin 9 n = 


kd 


equation (A.15) can be written 

£n(y) e Ti0«* _ ^ e ifc(v«n9„T*co»#«) _ g j <e( - V «in *« co* *„ ) J 


(A.16) 


(A.17) 


(A.18) 


because 


I. - 

- ^ w- 

= k\f 1 — sin 2 0 n 
= kcos9 n - 


(A. 19) 
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Figure A.2: Modal rays of the parallel plate waveguide. 


This form represents two crossing plane waves which make an angle of 0 n with the 
guide walls as shown in Figure A-2. 

The TM case is handled similarly, with the boundary condition being dH z /dy 
equal zero on the walls, 


H(x,y) 
(V 2 + k 2 )H(x,y) 


B(x,y) 

K(y) 

A. 

V X fi(x,y) 


E 



zH z (x,y) 

{h + ^ +kl ) H ‘ (x ' y) 


0 




zP n cos{—y) 



ju>eE(x,y ) 

k 


Z o 


Zp 

jk 



£ [A+e+(y)e + A n e n (y)e J0nX ] 


P n Zo 


, -A ,nir 
±V k cos (~d V ’ 


. m r . U7T 

x Wd’ m{ T v \ 


(A.20) 

(A.21) 
(A. 22) 

(A.23) 

(A. 24) 
(A. 25) 
(A. 26) 


(A.27) 
(A. 28) 
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(A. 29) 


/ 2 kYo 

V d(l n e n 

_ j 2 if n = 0 
_ 1 1 otherwi 


and the ray-optical form of equation (A.22) is given by 

h n (y)e** linX = z — -|- e jfc(-V*in»«T*co*<i.)j 


(A. 30) 


(A.31) 


Notice that for the TM case, the n = 0 term is included in the summation. This 
corresponds to a TEM mode. Ray-optically this mode is a plane wave traveling 
down the guide with $ n equal zero. 
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B. Orthonormal Modes of the Annular Waveguide 


This appendix derives the expressions for the orthonormal modal fields of the 
uniform annular waveguide. “Orthonormal” means that each mode carries power 
independently of all the other modes and they are normalized to carry unit power, 
as described more fully in Appendix A. The asymptotic ray-optical form of these 
modes is also derived. 

Figure B.l shows an annular waveguide of infinite extent. This makes no sense 
physically because the guide would join ends to form a ring. However, it is allowable 
mathmatically to assume that <f> goes to infinity in both directions. It is necessary 
to make this assumption to get the correct waveguide modes to describe the fields 
in small sections of an annular guide, such as the ones used to make up the inlet. 
Starting with the TE case, the electric field is in the z-direction and is a function 
of p and <j> 

E(p,<j>) = zE z {p,<j>). (B.l) 

The .E-Field is a solution to the wave equation 

(' V 2 + k 2 )E(p,<f > ) = 0 (B.2) 


which in scalar form is (using the cylindrical form of the Laplacian) 

E z {p, <f>) = 0 


1 d , 8 ^ Id 2 , 2 
+ * 


(B.3) 


pdp r dp' ' p'dtf 

where k is the free space wave number, 2ir/\. The boundary condition is that E z 
equals zero on the walls of the guide. Separating variables and substituting into 
equation (B.3), multiplying by p 2 and dividing out E z {p,<f>) gives 


dR(p) 


R(p) dp f dp J $(<£) d<f> 2 


+ 


l d 2 ${<f>) 


+ ( kp ) 2 = 0 


(B.4) 
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where the separation equation is 


eapa) = mm- (b-5) 

Introducing the separation constant v gives the two independent differential equa- 
tions 


i vm 3 

$(<£) d<f > 2 

(B.6) 

p dp \p dp ^ + KM 2 v 2 )R{p) = 0- 

(B.7) 

The solutions to equation (B.6) are 


*{4>) = e^ande'^. 

(B.8) 

(B.7) is Bessel’s equation and has solutions 


R(p) = Hl l \kp) and H*\kp) 

(B.9) 


where H^\kp) and H^\kp) are the Hankel functions of the first and second kinds, 
respectively, of order u and argument kp. Using linear combinations of these so- 
lutions and enforcing the boundary conditions gives the complete solution for the 
E-field 

EM) = f (B.10) 

n - 1 

where 

Mp) = zA n R n (kp) (B.ll) 

*n (kp) = l~. {H^kb^lkp) - (B.12) 

in which R n (kb) = 0 has been employed and where A n is the normalization con- 
stant of the n th mode and are arbitrary constants. If <f> had a periodic boundary 
condition corresponding to a 360 degree annular ring, the eigenvalues u n would be 
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completely determined. However, for the annular guide of infinite extent, the eigen- 
values must be found using the remaining radially p dependent boundary condition. 
Thus, they are found from the transcendental equation 


Rn(ka) = 0 (B.13) 

which usually must be solved numerically. 

As in a parallel plate guide, the summation generally only includes values of n 
for which u n is real. This is because for some values of n, u n will be purely imaginary 
and the mode will die out exponentially along the sods of the guide. 

The H -field of the guide is derived directly from the .E-field using one of Maxwell’s 
equations 

V x E(p,<f>) = -jupH( Pj <t>) (B.14) 

up = (B.15) 

Jo 

which gives 

SM) = 

= E [BfKip)'-*'* + fCMrte'- 4 ] (B.16) 

n= 1 

ht(p) = A„Y a ±^K(kf.) - ■ (B.17) 

A n is found by normalizing to unity the power carried by the n th mode. This 
done by integrating the Poynting vector over a cross-section of the guide and setting 
this equal to unity, 

= i. (B.i8) 

Substituting and changing the variable of integration to u = kp gives the normal- 
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ization coefficient 


j f^ a jj 2 t \ 

An ~ kzjkb RM u 




(B.19) 


The integral of equation (B.19) can be written in closed form using a property of 
Hankel functions. However, the resulting equation is quite long and cumbersome 
and contains derivatives with respect to the order u n of the Hankel functions. It 
was found that it is easier to just do the integration numerically. 

The ray-optical form of the waveguide modes in the annular guide is obtained 
using the large argument approximation to the Hankel functions (also known as the 
Debye Tangent approximation). For the Hankel function of the second kind, it is 
given by 


*?>(«) 




irx sin 7 6 


— i(* »in7— i/-y) 


where 


cos 7 


v 

x 


(B.20) 


(B.21) 


Using this approximation, the Hf£ term of a mode with the <f> variation included 
can be written as 




>2 


?rA:psin7 






cos 7 = - — . 

kp 


Changing variables gives the ray-optical form 

Hll ) {kp)e- jt ' n + « /^ e ~* fc “ +l ' n * 0) 


u = psmj 

< t > = 7 + ^0 


(B.22) 

(B.23) 

(B.24) 

(B.25) 

(B.26) 


where <f> 0 is a constant. This is the form of an outward traveling cylindrical wave. 
Figure B.2 shows the geometrical significance of this derivation. Notice that u is the 
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distance from a concentric caustic of radius v n jk and 7 is the angular displacement 
from this caustic. The term is handled the same way, 


//<■'(*) 


f — J2 

nx sin 7 


j(z sin 7— 1/7) 






__Z^gj(fcu+«'n<M 

irku 


(B.27) 

(B.28) 


which is the form of the corresponding inward traveling (convergent) cylindrical 
wave. 

Figure B.3 illustrates the modal ray form of the modes of the annular waveguide. 
It consists of crossing cylindrical waves which share a concentric caustic of radius 
u n /k. The angles the rays make with the walls of the guide are obtained easily from 
the ray-optic geometry, 


= cos-'(^) (B.29) 

= co» _I (~). (B.30) 


It happens that the modal ray caustic can lie inside the guide, as shown in 
figure B.4. This special case is referred to as a “Whispering Gallery” (WG) mode 
[]. The significance of this type of mode is that almost all of the power in the 
mode is confined between the modal caustic and the outer wall of the guide. The 
propagating portion of the fields never “see” the inner wall (although the modal 
caustic radius u n /k depends on it) because the fields are evanescent between the 
caustic and the inner wall, as shown in figure B.5. 

The TM case is handled similarly, with the boundary condition being dH x /dy 
equal zero on the walls, 


H{pA) 
(v 2 + k’)ii(p,t) 




\d_ 

pdp 


(p tp )+ 


1 d 2 
p 2 d<t> 2 



0 


h z M) 


(B.31) 


(B.32) 
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1 iM) = + 

n=l 


(B.33) 

^n(^) ^ zA n Rn(kp) 


(B.34) 

Rn(kp) = Yj [•®£ ) '(“) fl £ ) (M - fl£ y (»)fiS?(M] 

(B.35) 

R' n (ka) = 0 


(B.36) 

V x B(pA) = i^ f -E(pA) 


(B.37) 

k 

= Zo 


(B.38) 

= E + KK(p)* in ‘* 

n=l 


(B.39) 

ii(p) = a„z 0 Up^ukp) - hKM 


(B.40) 

< - |£f «•>?]''' 


(B.41) 


Notice that the eigenvalues u n for this case are found from the transcendental equa- 
tion (B.36). The ray optical form for the TM case is the same as the TE case with 
the appropriate change of constants. 
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C. Diffraction Coefficients 


The two types of diffraction effects considered in this report are from a wedge 
and from a discontinuity in curvature. The geometry for the wedge diffraction 
coefficient is shown in Figure C.l. The non-uniform GTD diffraction coefficient 
(Keller’s form) is given by [6] 

1 - 1 




sin " 

n 


n = 2 


n\/]2iHe [cos * — cos(^^) 

WA 


+ R,,h 


cos l - cos(*±£) 


7T 


(C.1) 

(C.2) 


where “s” designates the “soft” (TE) case and “/i” designates the “hard” (TM) case 
and R, h is the appropriate plane wave reflection coefficient for the surface. For a 
perfectly conducting surface, R, is -1 and Rh is 1. For an absorber coated surface, 
R ti h is given by equations (4.44) and (4.48) where is R t h . 

The geometry for the diffraction coefficient of a discontinuity in curvature is 
shown in Figure C.2. The non-uniform GTD diffraction coefficient is given by [7] 




&2 — &i 

ky/jZirk 


1 + cos(<ft + ft) £ 1 + cos(^ - tf>') 

(cos (f> + COS 4>') 3 s ’ h (co8 <j> + cos <f>') 3 


(C.3) 
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D. Symmetry Relations of Matrices 


Refering to Figure D.l, it is seen that some of the junction matrices are equal 
due to symmetry. These relations are as follows, 


1^33] 

= \Ru] 

(D.l) 

[S»l 

= [«3l] 

(0.2) 

l« 3 <] 

= (^3.) 

(0.3) 

P*..] 

= 1^.3] 

(D.4) 

PM 

= Pi.] 

(D.5) 

py 

= PM 

(0.6) 


There is a reversal of coordinate systems at the junction between sections 3 and 
4. Therefore, [T 45 ] ^ [T 43 ] and [T 64 ] ^ [T 34 ] contrary to the symmetry relations of 
(D.5) and (D.6) for the other two junctions. However, it can be shown that there is 
a simple relationship between the elements of the matrices [T 34 ] and [Tg 4 ]. For the 
TE case it is given by 


T 3 r = -(-i) n r 5 r 


and for the TM case by 


(D.T) 


T£" = (-1 ) n is n (D.8) 

where T™ n is the mn th element of matrix [T^]. 

As shown in section 3.2.1, the two transmission matrices at a junction are re- 
ciprocal, 


py = [r„] r 

[T32] = [T 2Z ] t 
[T4 a] = [T*) T 
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(D.9) 

(D.10) 

(D.ll) 




Figure D.l: Geometry of the S-shaped inlet. 


Pm] = P«] T . (D.12) 

With these simplifications, only six of the original sixteen junction matrices need 
to be evaluated explicitly. They are [Jli 2 ], [#21] > [#34] > [R43], Pnb and Pm] • If the 
inner radii of the curved sections are equal, i.e., r 2 = r 4 , then 

[•R34] = [^12] (D.13) 

[Rv] = [Rn\ ( D - 14 ) 

Pm] = [T12] (D. 15) 

leaving only [R12], P?2i], and p\ 2 ] to be evaluated. 
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E. Aperture Integration 


Figure E.l shows the relevant geometry for the Aperture Integration (AI). The 
aperture field is replaced by equivalent Kirchhoff electric and/or magnetic surface 
currents which radiate a far field pattern. 

For the TE case, the E-field is found from the electric and magnetic vector 
potentials F and A respectively, as [9] 

B(p,0) = -V x P(pJ) - jkZ 0 A(p,<>) + [V ■ A(p,l)j\ (E.l) 

where 

***> - < E - 2) 
- 7mIo fMe ^' dv - (E - 3) 

Zo is free space impedance and Fo = l/^o is free space admittance. The equivalent 
surface currents for the TE case are given by 


M(y) = E ap {y)x(-x) 

= 2E?(y) x (-x) 

= -yE?(y) (E.4) 

J{y) = -xxH at> (y) 

= -*X [iHj-w + $»;-(»)] 

= -sn; r (y). (E.5) 


Substituting (E.2), (E.3), (E.4) and (E.5) into (E.l), doing the vector operations 
and keeping only the terms of order 1/ y/p gives the final AI formulation for the TE 
case as 

E(p,«) = (E.6) 
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Figure E.l: Geometry for Aperture Integration 


The TM case is found analogously as follows: 

= V x A(pJ) - jkY a F(p,«) + [V • /'(/>.«)] (E.7) 

- -jjml M{y)e ' dy (E - 8) 

p-jkp r d ^ 

- «/„ (E - 9) 

/(y) = 

= — z x zH° p (y) 

= yffr(y) (E.10) 

M(y) = E ap (y) x (-*) 

= [*J5“ p (y) + y^“ p (y)] x (-*) 

= zE“ p (y) (E.ll) 

H(p,6) = Jq [cos0#“ p (y) — Y 0 E“ p (yj\ e~* kvain9 dy. (E.12) 
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